University  of 

Nebraska 

Lincoln 


December  2, 2000 

U.S.  Army  Research  Office 
4300  South  Miami  Boulevard 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 
ATTN:  AMSRL-RO-RI  (Sylvia  Hall) 

Dear  Sylvia, 

Enclosed  are  the  Final  Progress  Report  for  DAAG55-98-0462.  This  report  has  been 
granted  a  60  extension  to  submitted  before  December  15, 2000. 

If  you  have  questions  regarding  the  report,  please  contact  me  at  the  following 
address: 

Dr.  H.  Walter  Yao 
AMD 

915  Deguigne  St.  Room  270  D 
P.O.  Box  3453,  MS  143 
Sunnyvale,  CA  94088-3453 

Phone:  (408)  749-3544 

FAX:  (408)749-5585 
Email:  walter.yao@amd.com 

*  *  *  *  *  *  **  *  *  *  *  if  *******  if  ******  *** 

Best  Regards, 


Professor  of  Electrical  Engineering 


Department  of  Electrical  Engineering 
209N  Walter  Scott  Engineering  Center 
P.O.  Box  88051 1 
Lincoln,  NE  68588-0511 
Phone  (402)  472-3771 
FAX  (402  472-4732 


snc  QtJAtsrs  ;i-  ?rs;meD+ 

20010116  134 


REPORT  DOCUMENTATION  PAGE 


Form  A 
OMB  Ni 


Approved 
JO.  0704-0188 


Public  Reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comment  regarding  this  burden  estimates  or  any 
other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations 
and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188,) 
Washington,  DC  20503.  _ 


1.  AGENCY  USE  ONLY  ( Leave  Blank)  j2T  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 


December  2, 2000 


Final:  15  July  1998  -  14  July,  2000 


4.  TITLE  AND  SUBTITLE 


Optical  Properties  of  GaN  and  Other  m-Nitride  Semiconductor  Materials 
Studied  by  Variable  Angle  Spectroscopic  Ellipsometry 


5.  FUNDING  NUMBERS 


DAAG55-98-1-0462 


6.  AUTHOR(S) 
Huade  Walter  Yao 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Nebraska-Lincoln 
Research  Grants  and  Contracts  Office 
303  Administration  Building 
Lincoln,  NE  68588-0430 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.  S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

.  I  \  -  g’L- 


ii.  Supplementary  notes 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as 
an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  the  documentation. 


12  a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  12  b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


1.  ABSTRACT  (Maximum  200  words) 

In  this  2-year  project,  we  have  determined  anisotropic  optical  dielectric  functions  of  GaN,  AIN  and  Sapphire.  It  is  for 
the  first  time,  the  ordinary  and  extraordinary  optical  constants  of  the  important  Ill-Nitride  materials,  GaN,  AIN  and 
sapphire,  have  been  determined  by  variable  angle  spectroscopic  ellipsometry  (VASE),  in  the  energy  range  of  0.75  eV  to 
6.5  eV.  Our  results  indicate  that  anisotropic  optical)  properties  are  important  feature  of  the  GaN  and  other  hexagonal 
structured  m-nitride  materials.  The  extiabrdinary  optical  properties  of  those  materials,  which  were  omitted  before, 
influence  in  many  ways  the  optoelectronic  devices  properties  made  from  GaN,  AIN  etc. 

We  have  also  discovered  a  method  to  characterize  the  ordinary  and  extraordinary  dielectric  optical  functions  by  VASE 
under  c-plane  direction.  Our  method  indicates  that  the  ordinary  optical  dielectric  functions  (E  J_  <c>)  can  be  precisely 
determined  by  the  isotropic  mode  VASE  measurements  at  small  angles  of  incidence,  e.g.,  between  20  and  40  degrees, 
while  the  anisotropic  dielectric  functions  (E  I  I  <c>)  can  be  detected  most-sensitively  at  large  angles  of  incidence  (near  the 
pseudo-Brewster  angle). 


14.  SUBJECT  TERMS 

Optical  Properties;  GaN;  AIN;  Sapphire;  Anisotropic  optical  response;  Ordinary  and 
extraordinary  dielectric  functions;  Variable  angle  spectroscopic  ellipsometry  (VASE). 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OR  REPORT 

_ UNCLASSIFIED _ 

NSN  75404)1-280-5500 


18.  SECURITY  CLASSIFICATION 
ON  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF  ABSTRACT 


Standard  Form  298  (Rev.2-89) 
Prescribed  by  ANSI  Std.  239- 
298-102 


MASTER  COPY:  PLEASE  KEEP  THIS  "MEMORANDUM  OF  TRANSMITTAL"  BLANK  FOR  REPRODUCTION 
PURPOSES.  WHEN  REPORTS  ARE  GENERATED  UNDER  THE  ARO  SPONSORSHIP,  FORWARD  A  COMPLETED 
COPY  OF  THIS  FORM  WITH  EACH  REPORT  SHIPMENT  TO  THE  ARO.  THIS  WILL  ASSURE  PROPER 
IDENTIFICATION.  NOT  TO  BE  USED  FOR  INTERIM  PROGRESS  REPORTS:  SEE  PAGE  1  FOR  INTERIM  PROGRESS 
REPORT  INSTRUCTIONS. 


MEMORANDUM  OF  TRANSMITTAL 


U.S.  Army  Research  Office 
ATTN:  AMXRO-ICA  (Hall) 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


■  Reprint  (Orig  +  2  copies)  □  Technical  Report  (Orig  +  2  copies) 

■  Manuscript  (1  copy)  ■  Final  Progress  Report  (Orig  +  2  copies) 

O  Related  Materials,  Abstracts,  Theses  (1  copy) 
CONTRACT/GRANT  NUMBER:  DAAG55-98-1-0462 

REPORT  TITLE: 

Optical  Properties  of  GaN  and  Other  Ill-Nitride  Semiconductor  Materials  Studied  by  Variable  Angle  Spectroscopic 
Ellipsometry 

is  forwarded  for  your  information. 

SUBMITTED  FOR  PUBLICATION  TO  (applicable  only  if  report  is  manuscript): 


Sincerely, 


OPTICAL  PROPERTIES  OF  GaN  AND  OTHER  III-NITRIDE 
SEMICONDUCTOR  MATERIALS  STUDIED  BY  VARIABLE  ANGLE 
SPECTROSCOPIC  ELLIPSOMETRY 


FINAL  PREGRESS  REPORT 


Huade  Walter  Yao 


November  26, 2000 

U.S.  AMRY  RESEARCH  OFFICE 

DAAG55-98-1-0462 


University  of  Nebraska 
Department  of  Electrical  Engineering 


APPROVED  FOR  PUBLIC  RELEASE: 
DISTRIBUTION  UNLIMITED 


THE  VIEWS,  OPINIONS,  AND/OR  FINDINGS  CONTAINED  IN  THIS  REPORT  ARE 
THOSE  OF  THE  AUTHOR(S)  AND  SHOULD  NOT  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  ARMY  POSITION,  POLICY,OR  DICISION,  UNLESS  SO 
DESIGNATED  BY  OTHER  DOCUMENTATION. 


1. 


List  of  Appendixes 


Reprints  of  publication  bi  Material  Rewarch  Society  Symposium  Proceeding; 

Zavada  SPIE  P  ^  ^  H°Ve'  A’M‘  Wow<hak,  P.P.  Chow,  J.M 

Zavada,  SPIE  Proceeding,  3621,  73-84  (1999). 

2.  Reprints  of  publication  in  Journal  of  Applied  Physics: 

Yan,  J.  ApplPhys^1 85,' ^717-672^(1999)^  (a~All°i]  Sin8le  C^to/s  H.  Yao,  and  C.H. 
3-  Reprints  of  publication  in  Material  Research  Society  Symposium  Proceeding- 
H ,  c. 

Mat.  Res.  Soc.  Symp.  Proc.  591,  (2000).  '  '  '  ^0W/  J-  Han,  J.  M.  Zavada, 

4.  Reprints  of  publication  in  SPIE  Symposium  Proceeding: 

Proceedings  3938,  113-123  (2000).  '  ’  '  Rosamond' p-  p-  Chow,  J.  M.  Zavada,  SPIE 

5.  Reprints  of  publication  in  Journal  of  Applied  Physics: 

VaraMe  Angle  Spectroscopic  ElUpsomet’y  Ym  H^YaTT FS"  De,mmned  % 
Wowchah,  P.  P.  Chow, M.  Zavada,  h  ££  ^346^3^“  A  M' 

6.  Preprints  of  manuscript  submitting  to  Journal  of  Applied  Physics: 

Wowchak,  P.  P.  Chow, M.  Zavadafj.  Appl.  A'  M' 

7-  Copy  of  cover  page  of  the  Ph.D.  Dissertation: 

Related  Materials  Semiconductors  and 

2.  Statement  of  problems  studied 

optica,  prop^te^ 


\“h“SPhire  In  tWS  ^  the  most  important 

A.  Ordinary  and  extraordinary  optical  dielectric  functions  of  GaN 

devi7he  Wjde  tand  gap  semicond“ctor  GaN  is  an  important  material  for  light  emitting 
device  applications  in  the  green,  blue  and  UV  regions.  GaN  fiw  8  8 

substrates  usually  have  hexagonal  structure,  which  is  optically  anisotropi^The  ^Phlre 

optica,  dielectric  functions  of  GaN  have  been  studied  LensLly  s  nee  960  s  BuZl 

extraordinary  dielectric  functions  of  GaN,  on  the  other  hand,  are  not  well  dented 

oarized  reflectance  spectroscopy  is  a  straightforward  technique  to  measure  the 

extraordinary  response  of  GaN.  However  it  may  be  lack  of  accuracy  because  not  oNy 

roughnessc  “  “  a*‘'a*  ^  ^  “  *  lack  °f  on  surface 

study  the  o  d  WOTk'  TaWe  anf’e  SpeCtrosc°Pic  elliPsometry  (VASE)  has  been  employed  to 
?r  Mi*  , ,  17  extraordina‘y  optical  dielectric  response  of  hexagonal  gallium  nitride 

MBE  R1"  7  8r0W"  °n  C"P’ane  SaPPhire  substrates  (“-A1A)  by  molecular  beam  epitaxy 

MBE).  Room  temperature  VASE  measurements  were  made,  in  the  range  of  0.75  ,0  6  5  eV  a" 

the  angle  of  incidence  in  between  of  20  and  80  degree.  VASE  data  simulations  by  isotropic 

and  anisotropic  models  indicate  that  the  ordinary  optical  dielectric  functions  (E 1  <c>)  can  be 

precisely  determined  by  the  isotropic  mode  VASE  measurements  at  angles  of  incidence 

be  ween  20  and  40  degrees  while  the  anisotropic  dielectric  functions  (E  I  I  <c>7‘an7e 

detected  most-sens, lively  at  large  angles  of  incidence  (near  the  pseudo-Brewster  angle).  Thus 

standar^VASE  meeXtra°rdTry  dieIeCWc  functions  are  precisely  determined  by 

e  range;f  075  to  6-5  ev- The  surface  *«*  a>» 
Sidered  in  the  VASE  data  analysis  in  order  to  separate  it  from  the  anisotropic  effect  The 

VASE  data  is  ana  yzed  by  a  model  dielectric  function  based  on  the  GaN  critical  poTnt 

s  ructure,  which  allows  for  a  non-zero  extinction  coefficient  k  below  the  band  gap  The 

icknesses  of  these  GaN  films  are  accurately  determined  via  the  analysis  as  well.  ? 

B.  Anisotropic  optical  respone  of  AIN 

Variable  angle  spectroscopic  ellipsometry  (VASE)  and  micro  Raman  scattering  have 

high  anil  ^ t  t0  7  the  °ptlCal  anisotr°Py  and  optical  constants  of  AIN  films  grown  at 
high  and  low  temperatures  (HT  and  LT).  The  AIN  films  were  grown  by  metalorganic  vapor 

P  ase  epitaxy  (MOVPE)  and  molecular  beam  epitaxy  (MBE)  on  c-plane  sapphire  (a-Al20 ) 

ubstrates,  respectively.  Anisotropic  optical  phonon  spectra  of  AIN  have  been  measured 

the  potoon  oTth80  d  ‘he  b0PHCal  3X15  <C>  °f  A'N  18  dther  perP“di“la^  «  parallel  ,0 
r  £  2  a  °nZer°  off-diagonal  elements  Aps  and  A  of  Jones 

matrix  m  the  reflection  VASE  (RVASE)  measurements  indicate  that  the  <0  of  AIN  is  slightly 

hTaInT  S"face  ,"0rmal  dUe  t0  SUbStrate  miscut  The  ordinary  optical  constants  of  both 
HT  AIN  have  been  determined  spectroscopic  ellipsometry  at  small  angles  of  incidence  so  that 


•  the  extraordinary  response  is  greatly  reduced.  The  film  thickness  along  with  the  surface 

overlayer  was  determined  via  the  VASE  data  analysis  as  well.  '  " 


C.  Amsotropic  optical  properties  of  Sapphire  (a-Al203) 


, '/7?h!T  ^  many  attractive  Properties  that  make  them  widely  applied  in 

so  ‘d-state  device  fabncations.  More  recently,  due  to  the  dramatic  progress  in  the  GaN 

related  blue  LEDs  and  lasers  field,  sapphire  as  a  good  substrate  candidate  becomes  even 

more  important  than  the  past  since  its  crystal  structure  is  similar  to  GaN  and  allows  to  realize 
a  small  lattice  mismatch  growth. 

Anisotropic  dielectric  functions  of  sapphire  have  been  determined  by  both  reflection 
and  transmission  variable  angle  spectroscopic  ellipsometry.  The  measurements  were  made 
on  a-plane  (2110)  sapphire  substrates  in  the  energy  range  of  0.75  eV  to  6.5  eV  at  room 
temperature.  The  orientation  of  the  optic  axis  of  the  a-plane  sapphire  sample  was 
determined  by  polarized  Raman  scattering  based  on  which  the  Euler  angles  are  set  for  the 
l  mg  mo  e  .  wo  pairs  of  Cauchy  user-defined  functions  were  constructed  to  describe  the 
optical  constants  of  both  ordinary  nx(oi)  and  extraordinary  n„(oo)  rays,  respectively.  The 
resulted  optical  constants  from  the  generalized  ellipsometry  analysis  are  in  good  agreement 
with  the  reported  data.  The  Kramers-Kronig  (KK)  relation  between  the  real  part  and  the 
imaginary  part  have  been  carefully  checked  for  both  ordinary  and  extraordinary  rays  The 
perfect  fitting  between  the  calculated  and  experimental  el(a»  function  indicates  that  both 
unctions  of  ordinary  and  extraordinary  are  KK  consistent.  The  refractive  index  difference 
etween  t  e  ordinary  and  extraordinary  rays  is  close  to  a  constant  (+0.008),  which  is  mainly 
determined  by  the  off-diagonal  signal  of  the  transmission  ellipsometry  data.  The  extinctiol 
coefficients  are  zero  below  6  eV,  and  increase  rapidly  above  6  eV. 

3.  Summary  of  most  important  results 

There  are  two  most  important  results  from  this  two-year  studies: 

1.  Anisotropic  dielectric  optical  constants  of  GaN,  AIN  and  Sapphire 

m  J!  •!  f°r  the  fir,St  time'  ^  °rdinary  and  extraordinary  optical  constants  of  the  important 
Ill-Nitride  materials,  GaN,  AIN  and  sapphire,  have  been  determined  by  VASE.  Our  results 

indicate  that  anisotropic  optical  properties  are  important  feature  of  the  GaN  and  other 
hexagonal  structured  Ill-nitride  materials.  The  - extraordinary  optical  properties  of  those 

materials,  which  were  omitted  before,  influence  in  many  ways  the  optoelectronic  devices 
properties  made  from  GaN,  AIN  etc. 

LiconduSrmaSth0d  °f  CharaC,erize  °Ptical  <*  biaxial  structured 

It  has  been  well  known  a  difficult  task  to  precisely  determine  the  optical  properties  of 
uniaxial  structured  semiconductor  materials,  especially  when  the  c/a  ratio  is  small.  We  have 


"  fiSTr1  a“;°  Characteri“  ordinary  and  extraordinary  dielectric  optical 
ftmcttons  by  VASE  under  c-plane  direction.  Our  method  indicates  that  the  ordinary  optical 
te  ectric  functions  (E  X  <c>)  can  be  precisely  determined  by  the  isotropic  mode  VASE 
measurements  at  small  angles  of  incidence,  e.g„  between  20  and  40  degrees  while  th! 
anisotropic  dielectric  functions  (E  I  I  <c»  can  be  detected  most-sensitively  a,  large  angles  o 

incidence  (near  the  pseudo-Brewster  angle).  8  S 
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ABSTRACT 

Generalized  variable  angle  spectroscopic  ellipsometry  (VASE)  and  Raman  scattering  have  been  employed  to  study  the  optical 
anisotropy  of  GaN/Sapphire  structures.  The  GaN  films  were  grown  by  hydride  vapor  phase  epitaxy  (HVPE)  and  molecular 

beam  epitaxy  (MBE)  on  both  m-plane  and  c-plane  sapphire  ((X-AI2O3)  substrates,  respectively.  Anisotropic  optical  phonon 
structure  of  sapphire  have  been  measured,  based  on  which  the  optical  axis  of  sapphire  substrate  has  been  determined.  A  541 
cm'1  TO  phonon  of  GaN  grown  on  m-plane  sapphire  substrate  has  been  discovered  experimentally  which  is  due  the  coupling 
of  Aj  and  Ei  TOs.  Optical  axis  orientation  of  GaN  film  on  m-sapphire  has  been  fully  determined  by  the  anisotropic  angular 

dependence  of  the  coupled  TO  phonon.  Off-diagonal  elements  Apst  and  Aspt  of  transmission  VASE  (TVASE)  are  very 
sensitive  parameters  related  to  the  optical  anisotropy.  The  optical  axis  orientation  of  GaN  on  m-sapphire  has  also  been 
accurately  determined  by  TVASE  at  two  special  sample  positions.  The  optical  anisotropy  due  to  GaN  film  and  sapphire 
substrate  has  been  successfully  separated  at  90°  sample  position  allowing  to  study  the  optical  anisotropy  of  GaN  film  only. 

Keywords:  GaN,  sapphire,  Raman  scattering,  generalized  ellipsometry,  optical  anisotropy 

1.  INTRODUCTION 

The  wide  band  gap  semiconductor  GaN  and  related  materials,  with  their  excellent  thermal  conductivity,  large 
breakdown  field,  and  resistance  to  chemical  attack,  have  a  very  promising  application  potential  for  both  high  temperature 
electronic  devices  and  short  wave-length  optical  emitters.1  2  The  recent  development  of  high-brightness,  blue  and  green  light 
emitting  diodes  (LEDs),3'  room  temperature  pulsed4  and  continuous-wave  (CW)5  quantum  well  lasers  has  greatly  encourage 
researchers  to  continue  the  work  on  these  materials.  Most  of  the  GaN  films  were  grown  on  c-plane  sapphire  substrates,  and 
they  usually  have  wurtzite  crystal  structure  (a-GaN),  which  is  anisotropic  (uniaxial).  Several  valuable  optical  property 

studies  have  been  carried  out  on  the  GaN/c-sapphire  structure  by  both  Raman  and  ellipsometry  .  However,  there  are  not 
much  data  available  about  the  optical  anisotropy  of  GaN  possibly  due  to  the  shortage  of  GaN  samples  grown  on  sapphire 
substrates  other  than  c-plane  orientations.  For  GaN/c-sapphire  case,  it  is  difficult  to  study  the  optical  anisotropy  of  GaN  film 
since  the  polarization  of  most  of  the  probe  light-wave  is  perpendicular  to  the  optical  axis  <c>.  For  GaN  grown  on  m-plane 
sapphire,  on  the  other  hand,  its  optical  axis  <c>  is  certain  degree  away  from  the  surface  normal  as  pointing  out  by  Matsuoka 

et  a/8',  which  provides  the  possibilities  for  the  optical  anisotropy  study.  It  was  found  that_the  crystal  orientation  relation 
between  GaN  and  m-plane  sapphire  is  (0113)/  (OlTO)  as  the  interface  plane,  and  [0332]  /  [21 10]  as  the  in-plane  orientation 
preference.  Most  importantly,  it  was  concluded  that  GaN  grown  on  m-plane  sapphire  has  a  smoother  surface,  lower 
background  carrier  concentration  and  stronger  PL  intensity  which  indicate  that  GaN  grown  on  m-plane  is  superior  than  on  c- 

plane  sapphire  substrates..  In  this  work,  based  on  our  previous  optical  anisotropy  study  of  sapphire9 ,  a  GaN/m-sapphire 
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sample  grown  by  HVPE  has  been  studied  by  Raman  scattering  and  generalized  ellipsometry.  The  optical  axis  orientation  of 
GaN  film  grown  on  m-plane  sapphire  has  been  fully  determined  by  both  anisotropic  coupled  TO  phonon  and  the  off-diagonal 
elements  of  transmission  Jones  matrix. 

2.  Theory 

2.1  Raman  Scattering 

It  is  well  known  that  sapphire  and  GaN  are  all  optically  anisotropic  (uniaxial)  materials  due  to  their  rhombohedral  and 
wurtzite  crystal  structures.  Group  theoretical  analysis  shows  that  the  irreducible  representation  for  the  optical  modes  of 

sapphire  is  u 

r  =  Aig  +  2Aiu  +3A2g  +  5Eg  +  4Eu,  (1) 

and  the  irreducible  representation  for  the  acoustical  modes  is  A2U  +  Eu.  Since  the  unit  cell  has  center-of-inversion  symmetry, 
all  vibrations  that  are  Raman  allowed  are  infrared  forbidden  and  vice  versa.  More  specifically,  two  Alg  modes  and  five  Eg 
modes  are  Raman  active  only,  while  two  A2U  modes  and  four  Eu  modes  are  infrared  active  only.  The  Aju  and  A2g  vibrations 
are  neither  infrared  nor  Raman  active. 


The  irreducible  representation  ofWurtize  GaN  optical  phonons  is:11 

T  =  Ai  (z)  +  2B1  +  E^x,  y )  +  2E2 ,  (2) 

for  phonon  propagating  along  or  perpendicular  to  the  optical  axis  <c>.  Where  x,  y,  z  in  parentheses  repent  the  directions  of 
phonon  polarization.  The  Ai  and  E\  modes  are  both  Raman  and  infrared  active,  two  E2  modes  are  only  Raman  active,  and 
Bi  modes  are  both  Raman  and  IR  silent.  Using  selection  rules  of  the  anisotropic  phonon  structures  of  sapphire  and  GaN,  the 
optical  orientation  can  be  quickly  determined. 

2.2  Generalized  Ellipsometry 


The  variable  angle  spectroscopic  ellipsometry  is  designed  to  accurately  determine  the  values  of  two  standard  ellipsometry 
parameters  \| /  and  A  ,  which  are  related  to  the  complex  ratio  of  reflection  (or  transmission)  coefficients  for  light  polarized 

parallel  (p)  and  perpendicular  (s)  to  the  plane  of  incidence.12  For  isotropic  material  systems, 

P  =  — -  =  tan(\| /)elA.  (3) 

The  electric-field  reflection  coefficient  at  an  incident  angle  of  <|>  is  defined  as  rp  (rs)  for  p  (s)-  polarized  light.  They  are  the 
diagonal  elements  of  Jones  matrix, 

Wsrople=[o  <4) 

The  \{/  and  A  are  not  only  dependent  on  material  dielectric  responses,  also  on  the  surface  condition,  sample  structure,  and 
other  properties  such  as  the  optical  anisotropy. 

For  the  anisotropic  material  system,  the  non-diagonal  elements  of  Jones  matrix  are  not  necessary  to  be  zero.  In  the 
transmission  VASE  configuration, 

ow-fe  Z}  <5) 

By  using  the  same  approach  with  considerably  more  algebra  involved,  we  can  still  predict  the  \| f  and  A  values,  but  the 
Fourier  coefficients  related  to  \j /  and  A  become  more  complicated.  The  generalized  ellipsometric  parameters  are  defined  as 
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Aspt  =  7^-  =  tan  \|/spelAsP 

tSS 

where  the  Apst  and  Aspt  describe  how  much  amount  p-  or  s-polarized  light  becomes  s-  or  p-polarized  light  after  the 
transmission,  respectively.  If  the  optical  axis  is  strictly  either  perpendicular  (perfect  c-plane  situation)  or  parallel  to  the 
electric  field  of  the  incident  beam,  the  Jones  matrix  is  diagonal.  While  it  may  not  true  for  GaN  films  grown  on  non  c-plane 
sapphire  substrates,  such  as  m-plane  sapphire  substrates  since  the  optical  axis  of  both  GaN  and  sapphire  may  form  an  angle 
respect  to  the  incident  polarization  so  that  the  off-diagonal  elements  of  Jones  matrix  may  not  vanish.  The  off-diagonal 
elements  Apst  and  Aspt  are  dependent  on  sample  position,  angle  of  incidence,  which  can  be  used  to  determine  the  optical  axis 
orientation  of  crystals.  Generalized  ellipsometry  was  first  introduced  by  Azzam  and  Bashara.14  The  recent  developments 
makes  this  technique  more  complete  and  powerful.15  ’ 16  In  a  word,  the  Generalized  ellipsometry  is  a  technique  which  can  be 
used  to  determined  all  the  elements  of  Jones  matrix  of  arbitrarily  anisotropic  and  homogeneous  layered  systems  with 
nonscalar  dielectric  susceptibilities. 

3.  Experiments 

The  Raman  spectra  were  taken  at  room  temperature  with  a  SPEX  1877E  triple  spectrometer  equipped  with  a  liquid-nitrogen 
cooled  CCD  camera.  The  excitation  light  source  was  an  Ar+  laser  operating  at  488  nm  with  the  output  power  kept  at  150 
mW.  A  back  scattering  geometry  was  employed  for  all  the  Raman  measurements. 

The  GaN  films  used  in  this  work  were  grown  by  both  HVPE  and  MBE  on  m-plane  and  c-plane  sapphire  substrates.  The 
backside  of  the  substrate  was  polished  so  that  the  transmission  type  VASE  measurements  could  be  made.  The  TVASE 
optical  measurements  of  anisotropic  mode  were  performed  in  the  energy  range  of  0.75eV  to  4eV  with  a  0.03eV  increment  at 
room  temperature  since  the  sample  is  opaque  above  the  GaN  band  gap  3.4eV.  A  rotational  sample  stage  allowed  us  to  take 
the  data  from  different  sample  positions  marked  by  the  angle  between  the  optical  axis  <c>  and  the  x-axis  of  the  lab 
coordinates.  Multiple  angles  of  incidence  from  -30°  to  +30°  were  used  in  order  to  see  the  angular  dependence  of  the  optical 
anisotropy  of  the  GaN/m-sapphire  sample. 

4.  Results  and  discussion 
4.1.  Optical  phonons  of  sapphire  and  GaN 

Raman  scattering  is  a  nondestructive  technique  used  to  detect  the  lattice  vibration  modes  related  with  crystal  orientation  and 
symmetry.  A  series  of  Raman  scattering  measurements  have  been  carried  out  on  c-plane  and  a-plane  sapphire  samples.  The 
spectrum  of  Raman  scattering  on  a  c-plane  sapphire  at  back-scattering  geometry  Z(Y,Y  +  X)Zare  shown  in  Fig.  1(a).  A 
similar  Raman  phonon  spectrum  shown  in  Fig.  1(b)  was  obtained  from  an  a-plane  sapphire  with  its  optical  axis  perpendicular 
to  the  polarization  of  the  incident  laser  beam.  While  a  different  spectrum  shown  in  Fig.  1(c)  was  obtained  from  the  same  a- 
plane  sapphire  with  the  optical  axis  parallel  to  the  polarization.  Only  two  TOs  appeared  for  this  configuration,  and  the  645 

cm'1  (LO,  Aig)  is  a  new  phonon  which  is  forbidden  in  both  Fig  1(a)  and  Fig  1(b)  configurations.  This  indicates  that  the 
645  cm'1  (LO)  mode  is  allowed  or  forbidden  when  the  optical  axis  of  the  a-plane  sample  is  parallel  or  perpendicular  to  the 
polarization  of  the  Raman  probe,  respectively. 

Optical  phonons  of  GaN  films  grown  on  c-  and  m-plane  sapphire  substrates  are  shown  in  Fig.  2.  The  dash  line  refers  to  the 
GaN  on  c-plane  sapphire,  a  E2  and  a  Aj  LO  phonon  are  observed  at  569  cm'1  and  734  cm'1,  respectively.  They  are 
consistent  with  reported  values.17  The  solid  line  represents  GaN  grown  on  m-plane  sapphire  substrate.  Except  the  E2  and 
E\  TOs,  a  new  541  cm’1  appeared  in  the  spectrum,  which  was  never  reported  as  the  experimental  result.  Notice  that  this 
phonon  is  not  any  of  the  standard  zone  center  phonons  propagating  along  or  perpendicular  to  the  optical  axis  <c>.n 
According  to  Loudon’s  theory,18'  this  one  belongs  to  a  coupled  mode  of  two  TO  phonons  (A\  TO  and  E\  TO)  when  the 
angle  0  between  phonon  propagation  direction  <q>  and  the  optical  axis  <c>  is  in  between  of  0°  and  90°.  For  non-zero  0 
values,  the  coupled  modes  satisfy  the  followings: 

co^TOl)  =  co2(A1(TO))sin2  0  +m2(E1(TO))cos2  9 

co2(T02)  =  co2(E1(TO)) 


(9) 

(10) 
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Fig.  1.  Optical  phonons  of  sapphire,  (a)  c-plane;  (b)  a-plane,  <c>s  1  Ej;  (c)  a-plane,  <c>s  ||  E,. 

where  the  TO]  is  the  new  coupled  TO  phonon  with  a  strong  0  angular  dependence,  and  the  co(T02>  is  same  as  cofE]  TO). 
Applying  equation  (9)  to  the  measured  541  cm'1  coupled  TO  phonon  mode,  the  related  angle  0  is  about  58°. 

In  the  back  scattering  geometry,  the  phonon  wave  vector  <q>  is  perpendicular  to  the  sample  surface.  In  order  to  located  the 
orientation  of  <c>  of  GaN  (<c>g),  the  reference  direction  <c>  of  sapphire  (<c>s)  has  to  be  determined  first.  Based  on  the 
analysis  of  Fig.  1,  we  can  use  the  645  cm-1  LO  phonon  of  sapphire  to  find  the  <c>s  position.  As  shown  in  Fig.  3,  The  645 
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Fig.  2.  Optical  phonons  of  GaN  films.  Dash  line  is  for  GaN  on  c-plane  sapphire;  solid  line  is  for  GaN  on  m-plane  sapphire.  541  cm'1 
is  the  coupled  TO  phonon  of  Ai  and  Ei  TOs. 

cm'1  LO  of  sapphire  is  forbidden  when  <c>s  perpendicular  to  the  polarization  (dash  line),  and  it  is  allowed  as  <os  parallel  to 
the  polarization  (solid  line).  Therefore  a  reference  plane  can  be  formed  by  <q>  and  <c>s,  which  is  perpendicular  to  the 
sample  surface.  The  task  now  is  to  determine  where  the  <c>g  locates  at,  either  in  the  reference  plane  or  in  the  plane  which  is 
perpendicular  to  both  sample  surface  and  the  reference  plane,  or  in  some  intermediate  plane.  This  problem  can  be  solved  by 
tilting  the  sample  to  change  the  angle  0.  If  <Og  is  in  the  reference  plane,  by  rotating  the  sample  about  the  x-axis  which  is 
perpendicular  to  the  reference  plane  will  result  in  a  maximum  angle  0  change,  which  should  create  a  different  coupled  TO 
phonon  according  to  the  equation  (9).  But  the  coupled  TO  frequency  will  keep  unchanged  when  rotate  the  sample  about  the 
<c>s.  In  fact,  this  is  the  case  for  the  GaN/m-sapphire  sample  we  have,  as  shown  in  Fig.  4  (a).  The  solid  line  refers  to  the 

Raman  data  of  the  untilted  sample  case,  which  has  the  coupled  TO  at  541  cm"1.  The  other  dash  lines  refer  to  the  -20°,  -30°, 
20°,  and  30°  tilted  sample  cases,  which  show  that  the  coupled  TO  phonon  moves  either  to  higher  or  lower  wave  number 
indicating  the  angle  0  is  changed  correspondingly  with  the  sample  rotations.  The  coupled  TO  wave  number  keeps 

unchanged  while  rotates  the  sample  about  the  <c>s  in  our  experiment.  The  calculated  and  experimental  values  of  the  coupled 
TO  phonon  frequency  as  a  function  of  angle  0  is  plotted  in  Fig.  4  (b).  Notice  that  the  relation  between  0  and  a  is  untrue  for 

0  =  58°  ±a,  (11) 

where  a  is  the  tilting  angle.  In  stead, 

0  =  58°±arcsin[sin(a)/rig]  (12) 
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Fig.  3.  Optical  phonons  of  GaN  film  grown  on  m-plane  sapphire  substrate  at  two  sample  positions  respect  to  the  polarization.  Dash 
line  represents  <c>s  JL  Ej  case  where  645  cm*1  LO  of  sapphire  is  forbidden,  and  solid  line  represents  <c>s  ||  Ej  case,  where  the  645 
cm*1  LO  of  sapphire  is  allowed. 

is  true  when  the  sample  is  tilted  away  a  degree.  Where  ng  =  2.4  is  the  refractive  index  of  GaN  at  the  probe  laser  wavelength 
(488nm).  By  taking  account  the  interface  refraction  effect  on  the  <q>  direction,  the  experimental  data  is  in  good  agreement 
with  the  theory. 

4.2.  Generalized  eilipsometry  of  GaN/m-sapphire 

Optical  anisotropy  can  also  be  measured  by  generalized  eilipsometry,  especially  the  transmission  type  measurements  for 
transparent  materials.  The  TVASE  data  of  an  a-plane  sapphire  are  shown  in  Fig.  5  as  a  good  example.  The  data  were 
obtained  at  0°,  1°  and  60°  sample  positions  in  the  range  of  0.75  eV  to  5.8  eV.  The  angle  of  incidence  are  20°  and  30°, 

respectively.  It  can  be  seen  clearly  in  Fig.  5  (c)  that  the  off-diagonal  elements  Apst  and  Aspt  of  the  Jones  matrix  have  large  \j / 
values  indicating  the  existence  of  the  optical  anisotropy  in  sapphire  crystal.  The  magnitude  of  the  off-diagonal  signal  in  Fig. 
5  (b)  is  still  in  reasonable  size  even  for  only  about  1°  difference  between  <c>s  and  x-axis.  The  zero  Apst  and  Aspt  in  Fig.  5  (a) 
are  obtained  at  the  sample  position  of  <c>s  ||  x  (the  same  result  for  <c>  ±  x).  By  finding  the  zero  Apst  and  Aspt  position,  the 

optical  axis  of  sapphire  can  thus  be  fully  determined.  The  results  from  Fig.  5  agree  each  other  very  well  with  Raman  analysis 
in  Fig.  1,  but  with  more  accuracy  (<1°).  The  one  drawback  of  TVASE  method  for  determining  the  optical  axis  is  that 
TVASE  can  not  distinguish  the  0°  and  90°position  since  all  the  off-diagonal  elements  are  zero  at  these  two  positions. 

In  order  to  reveal  the  optical  anisotropy  of  GaN,  a  series  of  anisotropic  mode  TVASE  measurements  have  been  made  on  the 
GaN/m-sapphire  sample,  as  shown  in  Fig.  6.  The  selected  typical  data  here  are  taken  from  four  sample  positions:  0°,  1°,  90° 
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Fig  4.  (a)  Anisotropic  angular  dependence  of  the  coupled  TO  phonon  of  GaN  grown  on  m-plane  sapphire.  The  sample  was  rotated 
about  the  x-axis  which  is  in  the  sample  surface  plane  and  perpendicular  to  the  sapphire  optical  axis  <c>s, 

and  92°,  as  shown  in  Fig.  6.  The  angle  a  is  between  <c>s  and  x-axis.  The  spectra  are  cut  off  at  about  3.4  eV  since  no 
transmission  above  the  GaN  band  gap.  Multiple  angles  of  incidence  (0°,  10°,  20°  and  30°)  have  been  used  for  collecting 
strong  off-diagonal  signal.  At  0°  position  shown  in  Fig.  6  (a),  the  off-diagonal  elements  Apst  and  Aspt  are  all  zero.  A  clear 
conclusion  can  be  drawn  from  Fig.  6  (a)  is  that  the  optical  axis  of  GaN  <c>g  must  lay  in  the  plane  of  incidence  somewhere 
between  x-axis  and  z-axis,  since  <Og  could  be  along  neither  x-axis  nor  y-axis  direction  for  the  zero  Apst  and  Aspt,  and  we 
already  know  from  Raman  scattering  that  the  angle  between  <Og  and  <Os  is  larger  than  zero  and  less  than  90°.  Recall  the 
calculation  result  based  on  the  coupled  TO  phonon  mode  at  541  cm'1  and  equation  (9),  the  angle  0  between  <q>  and  <c>s  is 
about  58°.  We  now  also  know  that  the  so-called  reference  plane  formed  by  <q>  and  <c>s  defined  previously  is  the  plane  cf 
incidence.  Therefore,  The  orientation  of  <c>g  can  be  fully  determined,  which  is  in  the  reference  plane  (or  the  plane  cf 
incidence)  and  32°  away  from  <Os  (or  x-axis).  To  test  the  accuracy  of  this  measurement,  another  TVASE  data  was  taken  at 
1°  position  as  shown  in  Fig.  6  (b).  The  sizable  Apst  and  Aspt  were  detected  at  all  the  angles  of  incidence  indicating  the 
sensitivity  of  TVASE  is  high  enough.  At  a  larger  angle  position,  much  bigger  pairs  of  Apst  and  Aspt  can  be  obtained,  which 
is  a  good  evidence  of  the  optical  anisotropy  of  the  whole  structure.  However,  it  is  difficult  to  know  how  much  anisotropy 
comes  from  GaN  at  these  positions.  Because  of  the  special  orientation  relationship  between  <c>g  and  <c>s  for  GaN  on  m- 
plane  sapphire,  the  90°  position  can  serve  very  well  for  this  purpose  as  shown  in  Fig.  6  (c).  The  nonzero  Apst  and  Aspt 
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Fig.  4  (b)  A  comparison  of  the  coupled  TO  positions  from  the  theoretical  calculation  and  the  experimental  data. 

appeared  in  Fig.  6  (c)  are  purely  due  to  the  anisotropy  of  GaN  film,  since  the  contribution  from  sapphire  substrate  to  the  Apst 
and  Aspt  is  zero.  Data  in  the  Fig.  6  (c)  is  very  valuable  especially  for  GaN  dielectric  functions  study.  Only  two  degree  away 
from  the  90°  position,  the  Apst  and  Aspt  are  modified  by  the  sapphire  substrate  dramatically  as  shown  in  Fig.  6  (d).  This 
indicates  that  the  optical  axis  orientation  of  GaN  can  be  determined  accurately. 

The  angle  of  incidence  dependence  of  the  optical  anisotropy  is  showed  in  Fig.  7  (a)  and  (b)  refer  to  two  sample  positions,  0° 
and  90°,  respectively.  The  TVASE  data  was  taken  at  three  photon  energies  of  1.0  eV,  1.6  eV,  and  2.2  eV  (arbitrary  chosen 
below  GaN  band  gap).  The  angular  increment  is  1°.  In  Fig.  7  (a),  the  off-diagonal  elements  are  zero  in  the  whole  angular 

range.  This  is  due  to  that  the  <c>g  and  <c>s  are  all  in  the  plane  of  incidence.  In  Fig.  7  (b).  A  different  behavior  is  observed 
for  the  same  three  wavelengths  below  the  band  gap.  This  is  because  that  the  optical  axis  of  GaN  is  in  the  plane  perpendicular 
to  both  the  sample  surface  and  the  plane  of  incidence.  As  discuss  previously  the  magnitude  of  Apst  and  Aspt  are  all 
contributed  by  GaN.  The  Apst  and  Aspt  are  zero  at  normal  incidence  at  90°  position  which  is  again  the  proof  of  <c>g  in  the 
reference  plane.  The  bigger  Apst  and  Aspt  values  at  larger  angle  of  incidence  suggest  that  a  large  angle  of  incidence  is  usually 
a  better  choice  for  measuring  the  optical  anisotropy  by  TVASE. 

According  to  the  data  analysis  of  both  Raman  scattering  and  TVASE,  the  growth  configuration  of  this  GaN/m-sapphire 
sample  is  (0113) /  (0001)  as  the  interface  plane,  and  [0332] /  [01 10]  as  the  in-plane  orientationpreference,  which  is  different 
from  the  reported  structure  configuration  (0ll3)  /  (0110)  as  the  interface  plane,  and  [0332]  /  [21 10]  as  the  in-plane 
orientation  preference,  we  consider  this  is  the  other  option  for  GaN  grow  on  m-plane  sapphire  substrates. 
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Fig. 5.  Transmission  VASE  data  of  an  a-plane  sapphire  at  the  angle  of  incidence  of  20°  and  30°,  respectively,  (a)  oc-0°;  (b)  a-l°;  (c) 
a=60°.  a  is  the  angle  between  the  optical  axis  and  x-axis  of  the  lab  coordinates. 

5.  Summary 

Optical  anisotropy  nature  of  GaN  films  grown  on  m-plane  sapphire  substrates  was  characterized  by  Raman  scattering  and 
generalized  variable  angle  spectroscopic  ellipsometry.  A  541  cm'1  phonon  of  GaN  due  to  the  TO  coupling  was  observed 
experimentally  and  it  agrees  very  well  with  the  theoretical  prediction.  The  coupled  TO  phonon  of  GaN  film  grown  on  m- 
plane  sapphire  substrate  is  a  good  signature  of  its  optical  anisotropy,  and  the  angular  dependence  relation  can  be  used  to 
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Fig.  6.  Transmission  VASE  data  of  GaN/m-sapphire  at  the  angle  of  incidence  of  0°,  10°,  20°  and  30°,  respectively.  (a)<c>s  0°  offx- 
axis;  (b )<c>s  1°  off  x-axis;  (c)  <c>s  90°  off  x-axis;  (d)  <c>s  92°  off  x-axis.  Apst  and  Aspt  are  the  off-diagonal  elements  of  transmission 
Jonse  matrix. 


determine  the  location  of  <c>g.  It  is  also  found  that  the  off-diagonal  elements  of  transmission  type  Jones  matrix  are  very 
sensitive  to  both  the  sample  position  and  the  angle  of  incident.  The  optical  axis  of  GaN  is  fully  determined  by  using  the 
null  ApSt  and  Aspt  position  at  0°,  which  is  in  good  agreement  with  the  result  from  the  angular  dependence  Raman  scattering. 
According  to  the  analysis  of  both  Raman  and  TVASE,  the  <c>g  of  GaN  on  m-plane  sapphire  lies  in  the  reference  plane 
defined  by  <q>  and  <c>s,  and  it  is  32°  away  from  the  <c>s.  TVASE  measurement  at  90°  position  has  given  non-zero  Apst 

and  Aspt  which  are  purely  due  to  the  anisotropy  of  GaN  film,  this  kind  of  data  will  be  very  important  for  the  future  studies  of 
anisotropic  dielectric  functions  of  GaN.  A  new  growth  configuration  of  GaN  on  m-plane  sapphire  has  been  observed  as 
(0ll3)  /  (0001),  and  [0332]  /  [0110]  for  the  interface  plane  and  in-plane  orientation  preferences. 
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Anisotropic  dielectric  functions  of  sapphire  have  been  determined  by  both  reflection  and 
transmission  variable  angle  spectroscopic  ellipsometry.  The  measurements  were  made  on  cr-plane 
(2110)  sapphire  substrates  in  the  energy  range  of  0.75-6.5  eV  at  room  temperature.  The  orientation 
of  the  optic  axis  of  the  cr-plane  sapphire  sample  was  determined  by  polarized  Raman  scattering 
based  on  which  the  Euler  angles  are  set  for  the  fitting  model.  Two  pairs  of  Cauchy  user-defined 
functions  were  constructed  to  describe  the  optical  constants  of  both  ordinary  n±(a))  and 
extraordinary  nB( <*>)  rays,  respectively.  As  a  result  the  optical  constants  from  the  generalized 
ellipsometry  analysis  are  in  good  agreement  with  the  reported  data.  The  Kramers -Kronig  (KK) 
relation  between  the  real  and  the  imaginary  part  have  been  carefully  checked  for  both  ordinary  and 
extraordinary  rays.  The  perfect  fitting  between  the  calculated  and  experimental  6](co)  function 
indicates  that  both  functions  of  ordinary  and  extraordinary  rays  are  KK  consistent.  The  refractive 
index  difference  between  the  ordinary  and  extraordinary  rays  is  close  to  a  constant  (+0.008),  which 
is  determined  mainly  by  the  off-diagonal  signal  of  the  transmission  ellipsometiy  data.  The 
extinction  coefficients  are  zero  below  6  eV,  and  increase  rapidly  above  6  eV.  ©  1999  American 
Institute  of  Physics.  [S0021-8979(99)01409-7] 


I.  INTRODUCTION 

Sapphire  crystals  have  many  attractive  properties  that 
make  them  widely  applied  in  solid-state  device  fabrications. 
More  recently,  due  to  the  dramatic  progress  in  the  GaN  re¬ 
lated  blue  light  emitting  diodes  (LEDs)  and  lasers  field,  sap¬ 
phire  as  a  good  substrate  candidate  has  become  even  more 
important  than  in  the  past  since  its  crystal  structure  is  similar 
to  GaN  and  allows  us  to  realize  a  small  lattice  mismatch 
growth. 

It  is  essential  to  have  precise  measurements  of  the  opti¬ 
cal  constants  of  sapphire  for  further  optical  study  on  GaN 
and  related  materials,  since  the  optical  dielectric  response  is 
of  utmost  importance  for  optoelectronic  device  design.  Be¬ 
cause  of  the  rhombohedral  structure  of  both  GaN  and  sap¬ 
phire  crystal,  they  are  optically  anisotropic  (uniaxial).  To 
study  the  optical  properties  of  GaN  film  grown  on  sapphire 
substrate,  it  is  necessary  to  have  the  knowledge  of  the  optical 
response  of  sapphire.  This  is  the  main  purpose  of  our  study 
on  sapphire.  Although  the  optical  properties  of  sapphire  have 
been  studied  since  19581  and  continue  to  receive 
attention,2-4  most  of  the  results  are  about  the  ordinary  optical 
response  (the  optical  axis  is  perpendicular  to  the  light  polar¬ 
ization).  The  current  state  of  the  optical  measurements  of  the 
extraordinary  response  is  still  less  than  adequate;  only  sev¬ 
eral  discrete  values  are  available  which  were  obtained  by  the 
traditional  prism  method.5  Another  work6  does  present  the 
extraordinary  response  in  a  fairly  wide  energy  range,  but  the 
range  of  experimental  data  is  limited  to  1. 2-5.4  eV,  and  only 
the  normal  reflection  ellipsometry  was  employed  in  their 
work,  which  does  not  have  very  high  sensitivity  for  the  de- 
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termination  of  the  deference  of  the  refractive  index  between 
the  ordinary  and  extraordinary  rays.  Therefore,  more  precise 
optical  measurements  are  needed  to  determine  accurately  the 
anisotropic  response  in  a  wider  photon  energy  range. 

In  this  work,  we  report  the  determination  of  room- 
temperature  anisotropic  optical  constants  of  sapphire  crystal 
in  the  energy  range  of  0.75-6.5  eV  by  intensity  transmission, 
reflection  variable  angle  spectroscopic  ellipsometry  (R- 
VASE),  and  transmission  VASE  (T-VASE).  All  the  mea¬ 
surements  were  taken  from  a-plane  sapphire  samples.  The 
optical  constants  were  extracted  from  VASE  analysis  of  all 
three  types  of  measurements  by  a  multiple-sample,  multiple- 
model  method. 


II.  EXPERIMENTS 

Two  pieces  of  a-plane  (2lT0)  sapphire  substrates  were 
used  in  our  optical  measurements.  A  single-side  polished  one 
was  for  R-VASE,  and  the  other  double-side  polished  one 
was  for  intensity  transmission  and  T-VASE  measurements. 
For  these  a-plane  samples,  the  optic  axis  (c)  is  parallel  to  the 
sample  surface,  and  it  will  be  shown  later  how  the  (c)  orien¬ 
tation  is  determined  by  Raman  scattering. 

The  optical  measurements  were  performed  in  the  energy 
range  of  0.75-6.5  eV  with  a  0.03  eV  increment.  The  inten¬ 
sity  transmission  data  were  taken  under  both  (c)J_Ej  ((c)  is 
the  optic  axis,  and  Ej  is  the  electric  field  direction  of  the 
incident  light)  and  (c)HEj  conditions.  The  optic  axis  was  set 
at  45°  horizontal  for  the  R-VASE  measurements.  A  series  of 
angle  of  incidences  were  used  from  65°  to  75°.  The  most 
important  T-VASE  measurements  were  taken  at  three  sample 
positions,  the  angels  between  (c)  and  Ej  being  1°,  60°,  and 
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90°,  respectively.  Two  angles  (20°  and  30°)  of  incidence 
were  used  at  each  sample  position. 

Polarized  Raman  scattering  was  used  as  a  complemen¬ 
tary  tool  to  determine  the  optic  axis  orientation  of  the 
a-plane  sapphire  samples.  The  probe  light  was  a  488  nm 
argon  laser  beam.  The  scattered  light  was  detected  by  a 
charge  coupled  device  (CCD)  camera  under  a  backscattering 
geometry. 


III.  ELLIPSOMETRY  BASICS 


A.  Standard  ellipsometry 

The  variable  angle  spectroscopic  ellipsometry  is  de¬ 
signed  to  accurately  determine  the  values  of  two  standard 
ellipsometry  parameters  ip  and  A,  which  are  related  to  the 
complex  ratio  of  reflection  (or  transmission)  coefficients  for 
light  polarized  parallel  (p )  and  perpendicular  ( s )  to  the  plane 
of  incidence.7  For  the  isotropic  material  system 


Rp 

p=  — =  tan(^)g 


/  A 


(i) 


The  electric  field  reflection  coefficient  at  an  incident  angle  of 
<f>  is  defined  as  Rp(Rs)  for  p (s) -polarized  light.  They  are  the 
diagonal  elements  of  the  Jones  matrix  as 


[*^J  sample 


Rp 

0 


0 

Rs 


(2) 


The  ^  and  A  are  not  only  dependent  on  dielectric  functions, 
but  also  on  the  surface  condition,  sample  structure,  and  other 
properties  such  as  the  optical  anisotropy  of  the  material. 

The  pseudodielectric  function  (e)  can  be  obtained  di¬ 
rectly  from  the  measured  if/  and  A  values 


(e)  =  (e])  +  i(e2)  =  sm2  4> 


1  +  tan2  <p 


(3) 


In  the  case  of  an  air  ambient  over  a  bare  bulk  material  with 
an  ideal  smooth  surface,  the  pseudodielectric  and  intrinsic 
dielectric  functions  are  identical.  To  determine  optical  con¬ 
stants  of  a  sample  with  more  complicated  structure,  VASE 
data  must  be  analyzed  using  a  parametric  model  that  is  ad¬ 
justed  to  fit  the  measured  data.7 

For  the  anisotropic  material  system,  the  off-diagonal  el¬ 
ements  of  the  Jones  matrix  are  no  longer  zero.  In  the  reflec¬ 
tion  VASE  configuration,  it  can  be  written  as 


[^]  sample 


^pp  ^sp 
Rns.  R  s  c 


(4) 


By  using  the  same  approach  with  considerably  more  algebra 
involved,  we  can  still  predict  the  if/  and  A  values,  but  the 
Fourier  coefficients  related  to  if/  and  A  become  more 
complicated,8  and  similar  treatments  are  needed  for  the  trans¬ 
mission  VASE. 


B.  Generalized  ellipsometry 

In  principle  standard  ellipsometry,  which  is  related  to 
structures  which  reflect  or  transmit  p~  or  ^-polarized  light 
into  p-  or  s-polarized  light,  respectively,  can  be  used  to  de¬ 
termine  both  ordinary  and  extraordinary  constants  of  an  an¬ 


isotropic  sample  by  mounting  it  in  a  special  manner  so  that 
the  optic  axis  is  strictly  either  perpendicular  or  parallel  to  the 
electric  field  of  the  incident  beam.6,9  The  Jones  matrix  in  this 
case  is  strictly  diagonal.  But  there  are  several  problems  for 
this  application.  (1)  It  is  hard  to  make  a  perfect  alignment 
since  the  crystal  is  usually  cut  with  up  to  an  error  of  several 
degrees  so  that  the  off-diagonal  elements  of  the  Jones  matrix 
will  not  vanish  completely.  (2)  R-VASE  cannot  directly  give 
the  information  of  the  difference  of  n±  and  nN .  Especially 
for  a  very  small  n1-nlh  e.g.,  n0-ne= 0.008  for  sapphire,  it 
is  not  reliable  to  determine  the  n±  —  n^  only  by  the  standard 
R-VASE  data  which  are  dominant  by  the  diagonal  elements. 
(3)  For  most  material  analysis,  we  may  deal  with  thin  films 
grown  on  substrates,  like  GaN  on  sapphire.  This  is  the  sys¬ 
tem  which  has  an  anisotropic  material  on  another  anisotropic 
material,  and  even  if  we  know  the  orientation  of  the  optic 
axis  of  the  substrate,  the  optic  axis  of  the  GaN  film  is  prob¬ 
ably  unknown,  and  in  general  it  may  not  overlap  with  the 
substrate  optic  axis.  In  this  situation,  the  standard  R-VASE  is 
no  longer  adequate.  Therefore,  we  need  a  technique  which  is 
able  to  detect  the  off-diagonal  elements  and  also  has  a  cor¬ 
responding  model  to  analyze  the  data. 

Generalized  ellipsometry,  which  was  first  introduced  by 
Azzam  and  Bashara10  can  serve  this  purpose.  The  recent  de¬ 
velopments  make  this  technique  more  complete  and 
powerful.11,12  In  a  word,  the  generalized  ellipsometry  is  a 
technique  which  can  be  used  to  determine  all  the  elements  of 
the  Jones  matrix  of  arbitrarily  anisotropic  and  homogeneous 
layered  systems  with  nonscalar  dielectric  susceptibilities. 


C.  Data  type  and  modeling 


Within  the  generalized  ellipsometry,  several  new  terms 
have  been  defined.  For  reflection  ellipsometry,  the  elements 
of  the  Jones  matrix  are  AnE ,  Aps,  and  Asp 


[*^]sample 


AXpp  ivS] 

*ps 


=  R< 


11EE 

*ss 

^ps  ^pp 
L  ^pp  *SS 

AnE 

AnE 


'sp 


R * 


*sp 


(5) 


The  AnE  parameter  is  the  complex  ratio  of  reflection  coeffi¬ 
cients  similar  to  p  of  the  isotropic  material  case.  But  the  Aps 
and  Asp  are  the  ratios  counting  the  amount  of  /?-  (or  s -)  po¬ 
larized  light  that  converts  to  5-  (or  p-)  polarized  light  after 
the  reflection  from  an  anisotropic  material  system.  In  the 
generalized  VASE,  three  pairs  of  ij/  and  A  are  the  measured 
parameters,  which  associate  with  An£,  Aps,  and  Asp,  respec¬ 
tively.  In  the  same  manner,  the  Jones  matrix  elements  of 
transmission  VASE  are  defined  as  AnEn  Apsl,  and  Aspt,  re¬ 
spectively.  They  also  have  three  pairs  corresponding  to  if/ 
and  A  as  the  measured  parameters. 

In  order  to  obtain  the  optical  constants  by  analyzing  the 
experimental  data,  an  optical  model  containing  the  fitting 
parameters  that  can  fully  describe  sapphire  optical  property 
is  needed.  Since  sapphire  is  a  uniaxial  crystal,  i.e.,  it  has  an 
optical  axis  that  is  associated  with  the  optical  anisotropy.  In 
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FIG.  1.  The  Eulerian  angles  and  related  coordinates. 


addition,  since  sapphire  (a-Al203)  is  chemically  stable  in 
air,  no  surface  oxide  layer  was  considered  in  the  model. 
Therefore  an  one-layer  uniaxial  sapphire  model  can  be  con¬ 
structed  by  the  sample  thickness  and  two  sets  of  optical  con¬ 
stants,  ordinary  n±(a )),  kx(co)  and  extraordinary  n^w), 
^n(w).  Once  the  lab  coordinate  is  selected,  the  relative  spa¬ 
tial  orientation  of  the  optical  axis  of  the  sapphire  sample  can 
be  located  by  using  Euler  angles  <fi,  8,  and  if/,  as  shown  in 
Fig.  1.  Our  lab-fixed  axes  are  the  X0,  Y0,  Z0.  The  Z0  axis  is 
normal  to  the  sample  surface.  The  X0  and  Z0  axes  are  in  the 
plane  of  incidence.  The  Y0  axis  is  perpendicular  to  the  plane 
of  incidence.  The  sample  unit-cell-fixed  axes  are  X,  Y ,  and  Z. 
As  shown  in  Fig.  1,  Euler  angles  cf>,  8 ,  and  ^correspond  to 
rotations  about  the  Z  axis,  then  about  the  new  X  axis,  and 
then  again  about  the  new  Z  axis.  The  uniaxial  sapphire  model 
can  be  constructed  either  based  on  c-plane  or  0-plane  con¬ 
figurations.  If  an  a  plane  is  selected,  we  may  choose  the 
electric  field  of  the  extraordinary  ray  parallel  to  the  Y0  direc¬ 
tion  (or  X0  direction),  thus  the  X0  (Y0)  and  Z0  directions  are 
correspond  to  the  E  fields  of  the  ordinary  ray.  If  the  0-plane 
sample  is  mount  with  its  optic  axis  parallel  to  the  Y0  direc¬ 
tion,  all  the  Euler  angles  are  zero  under  such  conditions.  In 
general,  the  optic  axis  of  a  sample  may  lie  in  an  arbitrary 
direction  and  the  Euler  angles  have  nonzero  values.  For  in¬ 
stance,  if  0  =  20°  and  0=^=0°,  it  means  that  the  n±  of 
sapphire  is  orientated  20°  from  the  X0  axis  in  the  sample 
surface  plane,  the  nVi  is  orientated  20°  from  the  Y0  axis,  and 
the  nz~n±  is  still  normal  to  the  sample  surface. 

IV.  RESULTS  AND  DISCUSSION 

A.  Optic  axis  determination  by  polarized  Raman 
scattering 

In  our  study,  0-plane  substrates  (2 HO)  were  selected  in 
order  to  determine  the  optical  constants  of  both  ordinary 
and  extraordinary  n^o)).  As  mentioned  earlier  in  the 
modeling  part,  the  Euler  angles  in  the  model  depend  on  the 
location  of  the  optic  axis.  Therefore  how  to  determine  the 
orientation  of  the  optic  axis  becomes  very  important.  For  a 
given  small  piece  of  an  0-plane  sapphire  sample,  its  optic 
axis  can  be  determined  by  polarized  Raman  scattering.  The 
selection  rules  of  Raman  scattering  on  both  c-plane  and 
c-plane  samples  are  shown  in  Fig.  2  for  a  comparison.  The 
three  pair  spectra  refer  to  two  backscattering  geometry, 
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FIG.  2.  Polarized  Raman  effect  of  c-plane  and  c-plane  sapphire  samples:  (a) 
(c)lEit  c  plane;  (b)  (c)±Eh  a  plane;  (c)  <c)||Eh  a  plane. 


Z(Y ,Y)Z  (solid  line)  and  Z(Y,X)Z  (dashed  line).  Figure  2(a) 
is  from  a  c-plane  sample,  Figs.  2(b)  and  2(c)  are  from  an 
0-plane  sample,  but  with  a  90°  sample  rotation  between  Figs. 
2(b)  and  2(c).  The  similar  phonon  spectra  of  Figs.  2(a)  and 
2(b)  indicate  that  the  optic  axis  of  the  0-plane  sample  in  Fig. 
2(b)  geometry  is  perpendicular  to  the  electric  field  of  the 
incident  laser  beam  E,-.  The  new  strong  645  cm-1  longitu¬ 
dinal  optical  (LO)  phonon  peak  in  the  Fig.  2(c)  refers  that  the 
optic  axis  is  parallel  to  the  incident  E  field.  It  is  the  charac¬ 
teristic  phonon  peak  of  the  lattice  vibration  along  the  optic 
axis  direction.  Therefore,  by  either  minimizing  or  maximiz¬ 
ing  the  645  cm-1  peak  intensity  through  the  sample  rotation, 
the  optic  axis  of  an  0-plane  sample  can  be  easily  located  by 
polarized  Raman  scattering. 


B.  VASE  data  and  analysis 

Sapphire  is  a  colorless,  transparent  material  with  very 
low  reflectivity.  The  reflectivity  of  the  off-diagonal  elements 
are  too  small  to  detect.  The  reflection  VASE  data  of  both  a 
c-plane  and  an  0-plane  sample  are  shown  in  Fig.  3.  In  theory, 
both  Aps  and  Asp  should  be  zero  for  the  a  c-plane  sample,  but 
they  should  not  be  zero  for  an  0-plane  sample.  The  near  zero 
Aps  and  Asp  we  have  seen  in  Fig.  3(b)  is  simply  due  to  the 
weak  reflectivity  (about  1 X  10“6).  In  order  to  collect  all  the 
information  of  the  Jonse  matrix,  a  R-VASE,  three  T-VASE, 
and  two  intensity  transmission  measurements  have  been  per¬ 
formed  on  two  0-plane  sapphire  samples,  the  data  and  fit¬ 
tings  of  these  three  types  are  shown  in  Figs.  4,  5,  and  6, 


6720 


J.  Appl.  Phys.,  Vol.  85,  No.  9,  1  May  1999 


H.  Yao  and  C.  H.  Yan 


FIG.  3.  Reflection-VASE  data  of:  (a)  c-plane,  and  (b)  a-plane  ((c)  off  E, 
45°)  sapphire  samples. 


respectively.  It  can  be  seen  clearly  in  Fig.  5(b)  that  the  off- 
diagonal  elements  Apst  and  Aspt  of  the  T-VASE  data  have 
very  big  if;  values  indicating  the  existence  of  the  optical  an¬ 
isotropy  in  sapphire  material.  The  magnitude  of  the  off- 
diagonal  signal  as  shown  in  Fig.  5(a)  is  still  reasonable  in 
size  for  only  about  1  °  difference  between  (c)  and  E,- .  The 
zero  Apst  and  Aspt  in  Fig.  5(c)  are  obtained  at  the  sample 
position  of  (c)_LE,  (the  same  result  for  (c)IIE/).  The  optic 
axis  orientation  is  in  well  agreement  with  that  determined 
from  Raman  measurements.  The  nearly  identical  intensity 
transmission  spectra  of  Figs.  6(a)  and  6(b)  indicate  that  the 
difference  between  n±  and  nS|  is  very  small  even  when  using 
a  quite  thick  sample  (0.332  mm).  The  best  model  fit  results 
for  the  three  data  types  are  also  presented  in  the  solid  line  in 
Figs.  4,  5,  and  6,  respectively.  The  excellent  data  fit  indicates 


45°. 


FIG.  5.  Transmission  VASE  data  and  fittings  of  a-plane  sapphire:  (a)  (c)  off 
Ej  1°;  (b)  <c>  off  E,  60°;  (c)  <c)  off  E,  90°. 

that  our  model  correctly  describes  the  optical  anisotropy  of 
the  sapphire  material. 

In  the  uniaxial  sapphire  model,  the  known  parameters 
are  sample  thickness  which  was  measured  by  a  micrometer, 
and  the  Euler  angles  based  on  the  Raman  measurements.  The 


FIG.  6.  Intensity  transmission  data  and  fittings:  (a)  ordinary  ray  ((c)! E,) 
and  (b)  extraordinary  ray  ((c)  HE,). 
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TABLE  I.  Best  fitting  parameters  of  the  anisotropic  sapphire  optical  con¬ 
stants  in  the  energy  range  of  0.75-6.5  eV. 


A 

B 

C 

D 

E 

1.734  94 

0.025  60 

-0.008  24 

0.002  71 

-0.000  13 

*11 

1.727  00 

0.026  85 

-0.008  41 

0.002  71 

-0.000  13 

a 

b 

c 

d 

2.3279X  10~7 

0.281  23 

0.037  06 

0.000  00 

*1! 

4.3854X  10~7 

0.280  72 

0.043  70 

0.000  00 

four  unknown  quantities  are  n±  ,  k±  ,  and  ,  jfc1(  described  by 
two  user-defined  Cauchy  functions  which  are  shown  below: 

B  C  D  E 

n~A+i+>?+)?+\*'  (6) 

k  =  (\  X  10~5)X£Xexp{-[(\-£)/c]3}  +  d,  (7) 

where  X  is  wavelength  in  /im.  All  the  constants  obtained 
from  the  best  fitting  are  listed  in  Table  I.  Notice  that  we  did 
not  simply  make  the  k±  =  k{]  =  0  assumption  in  our  data  fit¬ 
ting  since  the  k  value  may  be  detectable  at  high  photon  en¬ 
ergy  close  to  the  band  edge  (about  9  eV)2,13  in  a  thick 
sample.  In  fact,  the  k  values  of  larger  than  1  X  10" 6  were 
measured  and  they  showed  a  non-negligible  effect  on  all  the 
transmission-type  data.  Therefore  the  k  =  0  assumption  may 
bring  a  large  error  on  the  n  results  above  6  eV,  since  the 
intensity  transmission  is  sensitive  to  sample  thickness  and 
nonzero  k  values.  With  an  accurate  thickness,  both  kL  and  k\\ 
can  be  accurately  determined  by  the  intensity  transmission 
spectra.  The  amplitude  of  n±  and  nl{  are  basicly  controlled  by 
both  AnE  and  AnEt,  and  the  difference  of  n±  and  nM  is  the 
key  source  of  the  big  Apst  and  Aspt  values. 


C.  Kramers-Kronig  consistency  check 


For  a  set  of  correct  complex  dielectric  functions  6= 

+  ie2  of  any  material,  its  real  and  imaginary  parts  by  nature 
obey  the  Kramers -Kronig  (KK)  relation:9,14 


Ci(E)-l 


2  r*E' 


E'e2(E’) 

7rdE'. 


(8) 


In  our  application,  due  to  the  limited  measurement  range  of 
0.75-6.5  eV,  a  modified  form  of  the  KK  relation  was  used  to 
check  the  KK  consistency  of  a  set  of  optical  constants: 


A, 


=  i  Ef-(EV 


~P  f 

ir  J  c 


6.5  eV  El  €2&QS{Ef) 
0.75  eV 


dE\ 


(9) 


For  a  specified  material,  the  KK  integral  is  numerically 
evaluated  to  calculate  €\  values  from  e2*  Two  zero- 
broadening  oscillators  and  a  fixed  offset  are  employed  to 
represent  the  spectral  response  outside  the  range  of  experi¬ 
mental  measurement.  Et  is  the  center  energy  location  of  a 
zero-broadening  oscillator;  and  A(  is  the  amplitude  of  the 
zero-broadening  oscillator  located  at  the  Et  position.  The  cal¬ 
culated  values  6^  were  compared  with  VASE  determined  by 
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FIG.  7.  Optical  constants  of  sapphire  for  both  ordinary  and  extraordinary 
rays  with  reference  data. 


6meas  trough  a  regression  analysis  by  varying  the  values  of 
A],  Ej,  and  until  calculated  and  measured  values 

match  as  close  as  possible. 

Figures  7(a)  and  7(b)  show  a  comparison  of  the  mea¬ 
sured  values  e™*™  and  the  calculated  values  ef  of  the  ordi¬ 
nary  ray  and  the  extraordinary  ray,  respectively.  Also  the 
fitting  parameters  of  KK  transformation  are  listed  in  Table  II. 
In  our  calculation  two  zero-broadening  oscillators  were  used 
outside  the  measured  region.  The  perfect  match  of  the  two 
curves  in  Fig.  7  indicates  that  the  dielectric  functions  of  sap¬ 
phire  obtained  from  our  VASE  analysis  are  KK  consistent. 

D.  Anisotropic  sapphire  optical  constants 

Two  pairs  of  sapphire  optical  constants  of  both  ordinary 
n±(a))  and  extraordinary  n,|(to)  rays  are  shown  in  Fig.  8.  It  is 
necessary  to  point  out  that  n±  is  slightly  larger  than  nl{ , 
which  is  consistent  with  other  reports.5,15  The  difference  be¬ 
tween  them  is  close  to  a  constant  value  0.008,  which  is  also 
in  well  agreement  with  References  5  and  15.  The  symbolized 
data  points  with  *  and  A  are  the  n±  and  values  from 
DeFranzo  and  Pazol,  in  which  the  traditional  prism  method 
was  employed.  It  also  can  be  seen  that  k  is  equal  to  zero  at 
the  energies  below  6  eV,  but  it  rapidly  increases  above  6  eV 
indicating  that  the  absorption  tail  can  start  well  below  the 
band  edge  (9  eV). 

V.  SUMMARY 

In  summary,  we  have  studied  the  anisotropic  optical  re¬ 
sponse  of  sapphire  single  crystal  by  generalized  ellipsometry 
with  three  different  types  of  measurement,  R-VASE, 
T-VASE,  and  intensity  transmission.  Polarized  Raman  scat- 


TABLE  II.  Fitting  parameters  in  the  Kramers -Kronig ’s  self-consistency 
check. 


Fitting 

parameters 

E,  (eV) 

ef1*1  (eV) 

Ordinary 

Pole  No.  1 

0.053  01 

0.462  27 

1.4061 

Pole  No.  2 

11.955 

31.363 

1.4061 

Extraordinary 

Pole  No.  1 

0.034  03 

0.734  74 

1.4299 

«ii 

Pole  No.  2 

11.929 

30.346 

1.4299 
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FIG.  8.  KK  consistency  check  of  the  dielectric  functions  of  sapphire:  (a)  KK 
fitting  of  the  ordinary  dielectric  function,  (b)  KK  fitting  of  the  extraordinary 
function. 

tering  as  a  convenient  and  nondestructive  tool  was  used  to 
determine  the  optic  axis  location  of  the  a-plane  sapphire 
samples.  For  the  transparent  sapphire  crystal,  the  T-VASE 
measurement  is  critically  important  for  collecting  entire  in¬ 
formation  of  the  Jones  matrix,  and  it  is  especially  helpful  for 
extracting  the  small  difference  between  n±  and  nu .  The  two 
sets  of  optical  constants,  both  ordinary  n±((o)  and  extraordi¬ 
nary  /i||(co),  in  the  energy  range  of  0.75-6.5  eV  resulted 
from  our  best  model  fit  are  in  well  agreement  with  the  avail¬ 
able  discrete  data  in  recent  reports.  The  KK  check  further 
confirms  that  the  resulting  dielectric  functions  of  sapphire 


are  Kramers-Kronig  self-consistent.  The  difference  of  the 
refractive  index  between  the  ordinary  and  the  extraordinary 
is  about  a  constant  of  0.008  in  the  range  of  0.75-6.5  eV.  The 
extinction  coefficients  of  both  ordinary  and  extraordinary 
rays  are  close  to  zero  below  6  eV,  but  they  increase  rapidly 
above  6  eV  showing  a  non-negligible  absorption  tail  below 
the  band  edge.  Combining  polarized  Raman  scattering  and 
the  generalized  ellipsometry  technique  could  allow  us  to 
study  the  optical  anisotropy  of  any  crystal  with  arbitrary  op¬ 
tic  axis  orientation. 
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ABSTRACT 

GaN  films  grown  on  GaAs  and  sapphire  substrates  by  molecular  beam  .epitaxy  (MBE)  and 
metalorganic  vapor  phase  epitaxy  (MOVPE)  at  both  low  and  high  temperatures  (LT  and  HT)  were 
characterized  by  Raman  scattering  and  variable  angle  spectroscopic  ellipsometry  (VASE).  Optical  phonon 
spectra  of  GaN  films  are  obtained  through  back-scattering  geometry.  Crystal  quality  of  these  films  was 
qualitatively  examined  using  phonon  line-width.  Phonon  spectra  showed  that  the  HT  GaN  has  wurtzite 
crystal  structure,  while  LT  GaN  and  GaN/GaAs  have  cubic-like  structures.  Thickness  nonuniformity  and 
defect-related  absorption  can  be  characterized  by  pseudo  dielectric  functions  directly.  Surface  roughness 
also  can  be  determined  by  using  an  effective-medium  approximation  (EMA)  over-layer  in  a  VASE 
analysis.  Anisotropic  optical  constants  of  GaN,  both  ordinary  and  extraordinary,  were  obtained  in  the 
spectral  range  of  0.75  to  6.5  eV  with  the  consideration  of  surface  roughness,  through  the  small  and  large 
angles  of  incidence,  respectively.  The  film  thickness  of  the  GaN  was  accurately  determined  via  the 
analysis  as  well. 

INTRODUCTION 

GaN  and  related  Ill-nitrides  wide  band  gap  semiconductors  have  attracted  much  attention  recently 
due  to  the  successes  of  blue  light  emitting  diode  (LED)  and  laser  diode  (LD)  fabrications  [1],  However, 
there  are  many  unsolved  problems  in  the  nitride  film  growth  field,  such  as  high  dislocation  density, 
cracks,  and  lack  of  surface  flatness  [2].  In  order  to  improve  the  film  quality,  more  advanced  and  sensitive 
film  characterizations  are  highly  required,  especially  in  nondestructive  ways.  By  measuring  the  phonon 
spectrum  of  a  film,  Raman  scattering  can  reveal  its  crystalline  quality,  crystal  structure,  alloy 
composition,  optical  anisotropy,  and  many  others.  Spectroscopic  ellipsometry,  on  the  other  hand,  not 
only  can  measure  the  film  optical  constants  and  its  thickness,  but  also  determine  the  surface  roughness, 
thickness  uniformity,  and  defect  related  absorption.  In  this  paper,  we  report  phonon  spectra  of  LT  and  HT 
GaN  films.  The  anisotropic  optical  constants  along  with  film  thickness,  surface  roughness,  and  thickness 
nonuniformity  of  the  HT  GaN  are  studied  by  VASE  as  well. 

THEORY 

It  is  well  known  that  a-GaN  is  optically  anisotropic  (uniaxial)  material  due  to  its  wurtzite  crystal 
structures.  Group  theoretical  analysis  shows  that  the  irreducible  representation  for  the  optical  modes  of 
Wurtize  GaN  [3] 

r  =  A1(z)  +  2B1  +  E1(x,y)  +  2E2,  (1) 

for  phonon  propagating  along  or  perpendicular  to  the  optical  axis  <c>.  Where  x,  y,  z  in  parentheses 
repent  the  directions  of  phonon  polarization.  The  Ai  and  Ei  modes  are  both  Raman  and  infrared  active, 
two  E2  modes  are  only  Raman  active,  and  Bi  modes  are  both  Raman  and  IR  silent. 

Spectroscopic  ellipsometry  is  another  powerful  nondestructive  method  that  can  be  used  to  study 
the  optical  and  structural  properties  of  crystals.  Two  standard  ellipsometric  parameters  t|/  and  A  are 
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related  to  the  complex  ratio  of  reflection  (or  transmission)  coefficients  for  light  polarized  parallel  (p)  and 
perpendicular  (s)  to  the  plane  of  incidence  [4].  For  isotropic  material  systems,  \|/  and  A  are  defined  as 

Rp  .. 

p  =  —  =tan(\|/)e  .  (2) 


The  electric-field  reflection  coefficient  at  an  incident  angle  of  <J>  is  defined  as  Rp  (Rs)  for  p  (s)- 
polarized  light.  and  A  are  not  only  dependent  on  dielectric  functions,  also  sensitive  to  the  surface 
conditions,  sample  structure,  and  other  properties  such  as  the  optical  anisotropy. 

The  pseudodielectric  function  <£>  can  be  obtained  directly  from  the  measured  \|t  and  A  values: 


<  e  >=<  >  +i  <  ej  >=  sin2  O 


1  +  tan  O 


1-P 
1  +  P. 


(3) 


For  the  case  of  an  air  ambient  over  a  bare  bulk  material  with  a  perfectly  smooth  surface,  the  pseudo 
dielectric  and  intrinsic  dielectric  function  are  identical.  To  determine  the  optical  constants  of  a  thin  film 
on  a  substrate,  VASE  data  must  be  analyzed  using  a  parametric  model  that  is  adjusted  to  fit  the  measured 
data.  A  regression  analysis  is  usually  used  to  vary  the  model  parameters  (e.g.,  optical  constants  or  layer 
thickness,  etc.)  until  the  calculated  and  measured  values  match  as  closely  as  possible. 


EXPERIMENTAL 


GaN  films  used  in  this  study  were  grown  on  c-plane  (0001)  sapphire  and  (001)  GaAs  substrates 
by  molecular  beam  epitaxy  (MBE)  at  high  (~700  °C)  temperature,  and  by  metalorganic  vapor  phase 
epitaxy  (MOVPE),  at  low  (~550  °C)  temperature.  These  samples  are  all  one-side  polished  to  meet  the 
need  of  reflection  ellipsometry  measurements. 

The  Raman  spectra  were  taken  at  room  temperature  with  a  SPEX  1877E  triple  spectrometer 
equipped  with  a  liquid-nitrogen  cooled  CCD  camera.  The  excitation  light  source  was  the  Coherent 
Innova  300  Ar+  laser  operating  at  514.5  nm  with  the  output  power  kept  at  300m W.  A  back  scattering 
geometry  was  employed  for  all  the  Raman  measurements.  Most  of  the  measurements  were  performed 
from  the  front  surfaces  with  the  focusing  spot  less  than  1  pm  in  diameter.  Micro-Raman  measurements 
were  also  made  on  the  cross-section  ,  so  that  the  anisotropic  optical  phonons  were  observed. 

The  standard  reflection  variable  angle  spectroscopic  ellipsometry  (RVASE)  measurements  were 
performed  in  the  spectral  range  of  0.75eV  to  6.5  eV  with  a  0.02eV  increment.  Small  angles  of  incidence 
(20°  and  40°)  were  used  to  determine  the  ordinary  optical  constants,  while  the  extraordinary  optical 
constants  of  GaN  were  determined  from  the  VASE  data  at  the  near-Brewster  angles  of  incidence  (60°, 
70°,  and  80°). 

RESULTS  AND  DISCUSSION 


GaN  films  grown  on  c-plane  (0001)  sapphire  substrates  usually  have  hexagonal  structure  and  its 
optical  axis  is  parallel  to  the  optical  axis  <c>  of  sapphire.  Optical  phonons  obtained  from  a  HT  GaN  film 
grown  on  c-sapphire  are  shown  in  Fig.  1  (a).  The  solid  line  represents  the  phonon  spectrum  taken  from 
the  sample  surface,  in  which  two  GaN  related  optical  phonons  (E2  569  cm'1  and  A1  LO  736  cm'1)  are 
observed.  In  this  case  the  optical  axis  of  GaN  film  is  perpendicular  to  the  laser  polarization.  The  418  cm’1 
and  750  cm'  phonons  are  due  to  sapphire  substrate.  When  the  laser  beam  was  focused  on  the  cross- 
section  of  the  GaN  film  and  let  the  optical  axis  of  GaN  (sample  surface  normal)  be  parallel  to  the  laser 
polarization,  another  set  of  phonons  was  obtained  as  indicated  by  the  dash  line  in  Fig.  1  (a).  The 
existence  of  A1  TO  (533  cm'1),  El  TO  (561  cm'1)  and  LO  (742  cm'1)  indicates  that  this  GaN  film  has 
hexagonal  structure.  The  narrow  line-widths  indicate  a  good  crystalline  quality.  Unlike  the  HT  GaN/c- 
sapphire  sample,  phonon  spectra  from  a  LT  GaN/c-sapphire  and  a  GaN/GaAs  show  that  they  all  have 
cubic-like  structures  [5],  as  shown  in  Fig.  1  (b).  The  broad  peaks  of  TO  554  cm'1  and  LO  732  cm'1 
indicate  that  the  crystal  defect  density  in  both  films  is  high  either  due  to  a  large  lattice  mismatch  or  the 
low  temperature  growth. 
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Fig.  1  Optical  phonon  spectra  of  GaN  films  obtained  under  backscattering  geometry  [  Z(Y, Y  +  X)Z  ].  (a) 
Phonon  spectra  from  a  hexagonal  GaN  film,  (b)  Phonon  spectra  from  cubic-like  GaN  films. 

Fig.  2  shows  the  pseudo  dielectric  functions  of  three  GaN  films,  a  1  pm  HT  GaN  on  c-sapphire 
(#1)  in  Fig.  2  (a),  a  2  pm  HT  GaN  on  c-sapphire  (#2)  in  Fig.  2  (b),  and  a  0.2  pm  LT  GaN  on  c-sapphire  in 
Fig.  2  (c).  The  solid  lines  represent  the  real  part  of  the  pseudo  dielectric  functions  of  each  film.  The  dash 
lines  represent  the  imaginary  part.  As  showed  in  Fig.  1  (a),  the  high  quality  of  the  #1  HT  GaN  film  can 
also  be  seen  from  its  pseudo  dielectric  functions  in  Fig.  2  (a).  Comparing  with  Fig.  2  (a),  the  #2  HT  GaN 
film  shows  some  different  properties.  The  severe  amplitude  damping  of  the  interference  fringes  on  both 
<£i>  and  <£?>  in  Fig.  2  (b)  indicates  the  poor  thickness  uniformity  of  this  film.  Also  the  average  value  of 
<£2>  in  the  near  below  band  gap  (3.4  eV)  region  is  much  larger  for  #2  HT  GaN  than  that  for  #1  HT  GaN. 
In  general,  the  <£2>  values  below  band  gap  for  a  perfect  GaN  should  be  zero.  This  high  <£2>  tail  may  be 
caused  by  the  crystal  defect  absorption.  An  even  higher  <£2>  tail  was  observed  in  the  LT  GaN  showed  in 
Fig-  2  (c).  This  is  not  a  surprise  for  a  low  temperature  grown  film.  The  non-damping  interference  fringes 
indicate  that  the  LT  GaN  film  has  good  thickness  uniformity. 

Surface  morphology  is  an  important  issue  for  III-V  semiconductor  epitaxial  layers,  especially  for 
Ill-nitrides.  Usually  the  thickness  of  surface  rough  layer  is  ranging  from  a  few  nm  to  tens  of  nm  [6,  7], 
Quick  and  nondestructive  measurement  of  the  surface  roughness  along  with  film  thickness  and  other 
optical  properties  is  the  big  advantage  of  spectroscopic  ellipsometry.  Ellipsometry  is  a  surface  sensitive 
technique  because  it  not  only  measures  the  light  intensity,  but  also  measures  the  phase  change  after  the 
reflection  (or  transmission).  To  study  the  incident  angle  dependence  of  surface  over  layer  thickness,  a 
series  of  model  simulations  were  performed  at  5.0  eV  (above  the  band  gap).  The  calculated  ¥  and  A 
spectra  are  showed  in  Fig.  3  (a)  and  (b).  A  three-phase  model  (surface  over-layer,  1.0  pm  GaN  film,  and 
c-sapphire  substrate)  was  used  in  the  simulations.  Anisotropic  optical  constants  of  sapphire  were  used  in 
the  calculations  [8],  Both  'F  and  A  spectra  at  near  the  pseudo  Brewster  angle  region  show  great 
sensitivity  to  the  roughness  layer  thickness.  With  the  increase  of  the  roughness 
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Fig.  2  Pseudo  dielectric  functions  of  GaN/c-sapphire  structures,  (a)  1  pm  high  quality  HT  GaN  on  c- 
sapphire;  (b)  2  pm  HT  GaN  on  c-sapphire;  (c)  0.2  pm  LT  GaN  on  c-sapphire. 

layer  thickness,  the  A  spectrum  as  a  function  of  incident  angle  no  longer  has  the  abrupt  turn  feature  near 
the  pseudo  Brewster  angle.  This  could  be  a  useful  sign  for  monitoring  the  surface  roughness.  The  full 
spectral  simulations  of  T  and  A  with  or  without  surface  roughness  layer  are  shown  in  Fig.  3  (c)  and  (d).  It 
shows  clearly  that  'F  and  A  are  much  more  sensitive  to  the  roughness  layer  thickness  in  the  above  gap 
region  than  below  gap.  It  is  due  to  the  small  penetration  depth  of  photons  with  energies  above  the  band 
gap  under  strong  absorption.  On  the  other  hand,  photons  with  energies  below  the  band  gap  can  almost  see 
the  whole  structure.  Therefore,  it  allows  us  to  measure  the  film  thickness. 
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Fig.  3  VASE  model  simulations  of  the  surface  roughness  effect.  dsur  is  the  surface  roughness  layer 
(EMA  with  50%  voids  and  50%  GaN)  thickness.  The  model  is  a  three-phase  structure  (surface  rough 
layer(dsur)/GaN  (l|im)/c-sapphire).  (a)  &  (b)  As  a  function  of  angle  of  incidence  at  5.0  eV;  (c)  &  (d)  As  a 
function  of  photon  energy  at  60°  and  80°  angles  of  incidence. 

The  ordinary  and  extraordinary  optical  constants  of  a  high  quality  HT  GaN  grown  on  c-sapphire  are 
shown  in  Fig.  4  (a).  The  results  were  extracted  from  a  multiple  angle  VASE  data  analysis  using  two 
parametric  semiconductor  models  for  GaN.  The  indices  of  refraction  of  extraordinary  are  about  3.5% 
higher  than  its  ordinary  counterpart  throughout  the  whole  spectrum.  The  zero  k  below  band  gap  indicates 
that  this  HT  GaN  film  has  very  low  defect  density.  A  comparison  of  optical  constants  of  a  HT  and  a  LT 
GaN  film  grown  on  sapphire  substrates  is  shown  in  Fig.  4  (b).  The  broadening  of  n  and  k  spectra  at 
critical  points  indicates  the  existence  of  high  densities  of  structural  defects  in  the  LT  GaN  film. 

SUMMARY 

In  summary,  we  have  studied  the  optical  properties  of  HT  and  LT  GaN  films  by  both  Raman 
scattering  and  variable  angle  spectroscopic  ellipsometry.  It  has  been  observed  that  HT  GaN  grown  on  c- 
sapphire  has  hexagonal  structure,  while  GaN  grown  on  GaAs  and  LT  GaN  on  c-sapphire  may  have  cubic- 
like  structure  that  strongly  depends  on  the  growth  conditions.  Defect  density  and  thickness  nonuniformity 
can  be  characterized  using  VASE  by  monitoring  the  pseudo  dielectric  function  line  shapes.  VASE  is  also 
a  surface-roughness  sensitive  technique  which  can  precisely  determined  the  surface  over  layer  at  near 
pseudo  Brewster  angle  region  in  the  above  band  gap  spectral  range.  The  ordinary  and  extraordinary 
optical  constants  of  a  high  quality  HT  GaN  were  fully  determined  by  a  multiple  angle  VASE  data 
analysis.  The  indices  of  refraction  for  extraordinary  ray  are  about  3.5%  larger  than  that  of  the  ordinary. 
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Fig.  4  Optical  constants  of  GaN  films  obtained  from  a  multiple  angle  VASE  data  analysis,  (a)  Ordinary 
and  extraordinary  optical  constants  of  a  HT  GaN  film;  (b)  Ordinary  optical  constants  of  a  HT  and  LT 
GaN  films. 
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ABSTRACT 

Variable  angle  specrroacopie  ellipsonaenv  (VASE)  and  micro  Ram.,  scncrn, 1  have  ^ 

<c>  of  AIK  is  either  perpendicular  o,  parallel  ro  rhe  polar, canon  of  the  inoidem  beam.  Konaero  off-d, agon,  elements  A* 

.Asp  of  Jones  mamiain"  e  rehecdon  VASE  (RVASE,  tT— — 

ZZZZZZ,  response  is  greatly  reduced.  The  Elm  thtchness  along  with 
the  surface  overlayer  was  determined  via  the  VASE  data  analysis  as  well. 

Keywords:  AIN,  sapphire,  Raman  scattering,  spectroscopic  ellipsometry,  optical  anisotropy,  optical  constants. 

1.  INTRODUCTION 

The  wide  band  gap  .AIN  along  with  GaN  and i  InN,  are  the  important 

devices  and  short  wave-length  optical  emitters.1-  AIN  films  grown  carried  out  on  textured5 

(a-AIN),  which  is  anisotropic  (uniaxial).  There  are  some  valuable  op  ica  p  p  rty  accurate  values 

and  single  crystal  AIN  films  by  various  kinds  of  techniques.  For  reflectance’can 

of  reflective  index  at  several  discrete  wavelengths  but  overl  since  it  does  not  contain  the  phase 

obtain  spectroscopic  data  but  do  not  have  enough  sensi  ty  _easuring  ^  spectroscopic  optical  constants  and 

information.  Ellipsometry,  on  the  other  hand,  is  an  ideal  choice  for  ^  g  P  .  £  wQrk)  ordinary 

structural  parameters  of  thin  films,  especially  when  the  films  ave  mt  r  transmission  variable  angle  spectroscopic 

optical  constants  of  both  HT  and  LT  ^  thin  turn  optical  reSponse  from  the  uniaxial  AIN  films  was 

ellipsometry  (TVASE)  and  reflection  VASE  (  )  .  ,  d  for  minimizing  the  error  caused  by 

greatly  reduced  by  using  small  angle  of  incidence.  Small  ang  e  o  inci  enc  interface  roughness  were 

the  substratemiscut  that  can  result  in  non-zero  off-diagonal  elements  of  Jones  matnx.  ^Surface  and  rate 
carefully  considered  in  the  VASE  model  analysis  using  the  effective  medium  approxima  ion 


"  Correspondence:  Email:  hyao@unl.edu;  Telephone:  402-472-5914;  Fax:  402-472-4732 


In  Light-Emitting  Diodes:  Research,  Manufartunng,  and 
E.  Fred  Schubert,  Editors,  Proceedings  of  SPIE  Vol.  3938  (auuu) 


2.  Theory 


2.1  Raman  Scattering 


The  irreducible  representation  of  Wurtize  AIN  optical  phonons  is:7 


of  phonon  polarization.  A,  and  E,  modes  are  polar  mo  s  ^  ^  ^  ^  Raman  and  infiared 

"  ll  - ~ —  -  -  —  —  *»  -  -  - 

active  Raman  phonon  modes  in  the  AIN  crystal. 

2.2  Generalized  Ellipsometry 

parallel  (p)  and  perpendicular  (s)  to  the  plane  of  incidence.  For  isotropic  material  systems, 

p  =  ^  =  tan(y/)e,A.  (-) 

The  electric^  reflection  coefficient  at  an  incident  angle  of*  is  defined  as  rp  (rs)  for  p  (s)-  polarized  light.  They  are  the 
diagonal  elements  of  Jones  matrix. 

The  *  »nl  not^  only  dependent  on  dielectric  fate, it*,  .bo  o.  the  soribce  condition,  s.ntpl.  and  othec 

properties  such  as  the  optical  anisotropy. 

susceptibilities. 

Fo,  the  anisotropic  tnateri.l  system,  the  „on-d, agonal  elements  of  Jones  matrix  ate  no,  necessaty  to  be  zero.  In  the  r=«=cnon 
VASE  configuration, 


1 sample 


rw  Tsv 

rps  rss 


By  using  the  same  approach  with  considerably  more  algebra  involved L,  w  e :  ca^n  Jp^eter^t^defined  as 

Fourier  coefficients  related  to  y  and  A  become  more  complicated.  The  generalized  ell.psom  P 

below. 
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Asp  =  JE-  =  taiuf/speAsP 
hs 

where  the  Aps  and  Asp  describe  how  much  amount  p-  or  s-polarized  light  becomes  s-  or  p-polarized  light  after  the  reflection, 
respectively.  If  the  optical  axis  is  strictly  perpendicular  (perfect  c-plane  situation)  to  the  electric  field  of  the  incident  beam,  the 
Jones  matrix  is  stricctly  diagonal.  While  it  may  not  true  for  real  AIN  films  grown  on  nominal  c-plane  sapphire  substrates 
since  the  substrates  are  usually  miscut  within  1°.  The  off-diagonal  elements  Aps  and  Asp  are  dependent  on  sample  positions, 
angles  of  incidence. 


3.  Experiments 

The  Raman  spectra  were  taken  at  room  temperature  with  a  SPEX  1877E  triple  spectrometer  equipped  with  a  liquid-nitrogen 
cooled  CCD  camera.  The  excitation  light  source  was  the  Coherent  Innova  300  Ar+  laser  operating  at  488  nm  with  the 
output  power  kept  at  150  mW.  A  back  scattering  geometry  was  employed  for  all  the  Raman  measurements. 

Two  identical  thin  AIN  films  (about  130  nm  thick)  were  grown  at  high  temperature  (1070°C)  side  by  side  by  MOVPE  on 
two  c-plane  sapphire  substrates  with  polished  and  unpolished  backsides  to  serve  the  transmission  type  and  reflection  type 
VASE  measurements.  About  20  nm  thick  LT  (525°C)  GaN  nucleation  layers  were  deposited  before .the  HT  AIN  growth. 
Another  two  films  were  grown  by  MBE  at  high  (800°C)  and  low’  (400°C)  temperatures  with  different  thickness,  respectively. 
The  very  thin  (-36  nm)  LT  AIN  film  was  grown  on  a  one-side  polished  c-plane  sapphire.  The  thickness  of  the  MBE  grown 
HT  AIN  film  is  about  1pm.  Both  isotropic  and  anisotropic  mode  VASE  measurements  were  performed  in  the  energy  ranee 
of  0.75eV  to  6.5  eV  with  a  0.02  eV  increment  at  room  temperature.  The  range  of  angle  of  incidence  is  from  20°  to  80°  for 
reflection  VASE,  and  0°  to  30c  for  transmission  VASE. 

4.  Results  and  discussion 

4.1.  Optical  phonons  of  AIN  films 

Raman  scattering  is  a  nondestructive  technique  used  to  detect  the  lattice  vibration  modes  related  with  crystal  orientation  and 
symmetry.  A  series  of  Raman  scattering  measurements  have  been  carried  out  on  various  kinds  of  AIN  samples.  The  optical 
phonon  spectra  of  a  1pm  thick  AIN  film  growm  on  c-plane  sapphire  at  back-scattering  geometry  Z(Y,Y  +  X)Z  are  shown 
in  Fig.  1  (a).  The  solid  line  w-as  obtained  from  the  front  surface,  in  w'hich  the  optical  axis  <0  of  AIN  (and  sapphire  substrate) 
is  perpendicular  to  the  incident  polarization  (Y).  Three  AIN  related  optical  phonon  were  observed  under  such  condition 
including  656  cm  3  E2  (high),  888  cm  1  Ai  (LO),  and  a  246  cm’1  E:  (lowf)  w'hich  is  not  shown  in  this  figure.  The  other  two 
phonons  of  577  cm  1  Eg  (LO)  and  750  cm’1  Eg  (LO)  are  due  to  the  sapphire  substrate.  The  dash  line  in  Fig.  1(a)  w-as 
obtained  from  the  side  cross-section  surface  in  w'hich  the  AIN  optical  axis  <c>  is  parallel  to  the  polarization  of  the  incident 
laser  beam.  There  are  three  new' phonons  appeared,  612  cm* !  Al  (TO),  670  cm*1  El  (TO),  and  909  cm*3  El  (LO).  All  the 
above  six  optical  phonons  measured  from  AIN  film  are  in  good  agreement  with  the  theoretical  calculation  based  on  the  rigid- 
ion  model.13  The  comparison  is  listed  in  Table  1.  Fig.  1  (b)  is  a  phonon  line-wddth  comparison  of  HT  and  LT  AIN  films. 
The  solid  line  represents  the  thin  LT  AIN  (36  nm)  grown  by  MBE;  the  dash  line  represents  the  thin  HT  AIN  (130  nm) 
growm  by  MOVPE.  The  broadening  of  AIN  E2  (high)  may  be  due  to  the  high-density  structural  defects  (such  as  group  HI 
vacancies)  in  LT  material  and  the  misfit  dislocations  due  to  the  large  lattice  mismatch.  It  is  observable  that  the  AIN  E2 
(high)  in  both  thinner  films  moves  to  lower  frequency  comparing  with  the  AIN  E2  (high)  in  Fig.  1  (a).  This  may  be  caused 
by  the  partial  strain  relaxation  in  stead  of  full  relaxation  in  a  thin  lattice-mismatched  film  even  for  a  thickness  of  130  nm. 


Table  1.  Zone-center  optical  phonons  in  wurtzite  AIN,  in  cm'1. 


Mode 

A,  (LO) 

A,  (TO) 

E,  (LO) 

E,  (TO) 

E2  (low) 

E2  (high) 

Theory12 

885  • 

614 

923 

668 

252 

660 

This  work 

888 

612 

909 

670 

246 

656 

Fig.  1.  Optical  phonons  of  HT  and  LT  AIN  films,  (a)  1pm  thick  MBE  grown  HT  AIN  film i  on  ^^'andLT  thin 

and  dash  line  represent  that  <0  is  perpendicular  and  parallel  to  the  incident  P°'  ^  A1^  film- dash  line  represents  130 

MOVPE  grown  AIN  films  on  c-plane  sapphire  substrates.  Solid  line  represents  36  nm  thick  LI  AIN  mm,  uasn  P 

nm  thick  HT  AIN  film. 


4.2.  VASE  data  analysis 


For  ideal  c-plane  samples,  its  off-diagonal  element  of  Jones  matrix  must  be  zero  as  we  pointed  out  in  the  theory  part.  But  for 
real  nominal  c-plane  samples,  since  there  is  a  small  miscut  (usually  within  1°)  on  the  substrate,  the  off-diagonal  elements 


Angle  of  incidence  (degree) 

Fig.  3.  Anisotropic  mode  RVASE  data  of  HT  AlN/c- sapphire  grown  by  MOVPE.  (a)  ¥  as  a  function  of  photon  energy  at  60° 
angle  of  incidence;  (b)  as  a  function  of  incident  angle  at  2.40  and  5.48  eV. 
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U  ■  P'o  ?  Fia  ■>  (a)  is  an  anisotropic  mode  RVASE  data  obtained  from  a  loO  nm 
may  not  vanish  completely,  as  shown  m  Fu.  2  =•  -  (  )  which  1S  near  the  pseudo  Brewster  angle.  It 

MOVPE  grown  AlN/c-sapphire  sample.  It  was  taken  at  60  an  nealisible.  The  incident  angular  dependent 

can  be  seen  clearly  that  the  off-diagonal  elements  ps  “ ‘  J"  shown  1  (b).  The  angular  increment  was  1°. 

measurements  were  made  at  certain  wavelengths  (L4 and  ^  ^  ^  ?hQC  they  „  negligible  at  the  angle  of  incidence 

The  off-diagonal  elements  are  detectable  m  e  an=  =  cross-conversion  of  p-  to  s-  or 

below  40°  or  above  80°.  Therefore,  small  angle  of  incidence  (<40  )  can  be 
s-  to  p-  polarization  caused  by  the  <c>  axis  offset. 

,  l ...  aic0  can  oreatlv  reduce  the  anisotropic  effect 

Small  angle  of  incidence  is  not  only  good  for  of  incidence,  the  ordinary  dielectric 

on  the  diagonal  elements  as  shown  in  our  previous  a  „  ^  element  without  losing' the  accuracy.  Thus 

Fig.  3  shows  the  nominal  sample  structures  and  the  by  MOWE.  The  nominal 

two  130  nm  AlN/c-sapphire  strucmres  on  one-side  and  double-side  polish^  Fio_  3  (b).is  the  VASE  model 

3* (al.n At' the  small  incident  angle  conduion.  all  the  layers  in  the  structure  are  simulated  as 
representing  the  struciuiw  m  i  *=•  ~ 


HT  AIN  130nm 


EMA  (AlN  +  37%  Void)  18  nm 

HT  AIN  118.4nm 


LT  GaN  20nm 


c-plane  Sapphire 


LT  GaN  6.5nm 

Sapphire 
035miii  i 

(b) 


LT  AIN  36nm 


LT  AIN  34.2nm 


10  mins  Nitridization 

c-plane  Sapphire 


EMA  (AlN+9.3%  AEO3)  22nm 

Sapphire 

0.35mm 

(d) 


j.  _  VAcp  model  (a)  Sample  structure  of  130  nm  thick  HT  AlN/c- 

grown  by  MBE.  (d)  VASE  model  of  the  LT  AlN/c -sapphire. 
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isotropic  material.  The  surface  overlayer  was  modeled  by  an  EMA  layer  to  account  for  the  surface  oxide  and  roughness.  The 
optical  constants  of  LT  GaN  and  sapphire  used  in  the  analysis  are  from  our  previous  study.  Fig.  3  (c)  is  the  nominal 
structure  for  36  nm  LT  AIN  grown  by  MBE.  Notice  that  there  is  a  10  minutes  nitrogenation  process  before  grow  the  LT 
AIN.  and  it  was  represents  by  a  EMA  layer  as  shown  in  Fig.  3  (d).  The  optical  constants  of  HT  and  LT  AIN  in  the  model 
were  represented  by  two  parametric  semiconductor  models,  which  are  constructed  based  on  their  critical  point  structures.  The 
details  about  the  parametric  model  is  described  in  reference.  15.  The  thickness  and  EMA  composition  listed  in  the  model  are 
the  best  values  from  VASE  data  fitting. 

The  VASE  data  and  best  fit  of  HT  AIN/c-sapphire  are  shown  in  Fig.  4.  Fig.  4  (a)  and  (b)  are  \|/  and  A  of  TVASE  at  four 
angles  of  incidence  (0°,  10°,  20°,  and  30°).  Fig.  4  (c)  and  (d)  are  the  \j/  and  A  of  RVASE  at  three  small  angles  of  incidence 
(20°,  30°,  and  40°).  The  two  sets  of  ellipsometric  data  were  coupled  together  during  the  data  fitting  to  further  reduce  the 
correlation  between  the  fitting  parameters.  The  resulted  AIN  layer  thickness  is  1 18.4  nm  with  18  nm  surface  overlayer.  The 
thickness  of  overlayer  is  consistent  with  previous  thermal  stability  study.16  Fig.  5  (a)  and  (b)  show  the  RVASE  data  and 


Photon  Energy  (eV)  Photon  Energy  (eV) 

Fig.  4.  VASE  data  and  best  fit  of  two  HT  AIN/c-sapphire  samples,  (a)  \| f  data  and  best  fit  of  TVASE;  (b)  A  data  and  best  fit 
of  TVASE;  (c)  y  data  and  best  fit  of  RVASE;  (d)  A  data  and  best  fit  of  RVASE.  Solid  lines  represent  experimental  data. 
Dash  lines  represent  the  best  model  fit. 
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20 


Fig.  5.  RVASE  data  and  best  fit  of  LT  AIN/c-sapphire.  (a)  vy  data  and  best  fit;  (b)  A  data  and  best  fit.  Solid  lines  represent 
experimental  data.  Dash  lines  represent  best  model  fit. 

best  fit  of  the  LT  AIN/c-sapphire.  The  resulted  film  thickness  of  LT  AIN  is  34.2  nm  which  is  very  close  to  the  nominal 
value  36  nm.  A  22  nm  interfacial  layer  was  formed  due  to  the  10  minutes  nitrogenation.  No  surface  overlayer  is  needed  in 
the  LT  AIN  VASE  data  analysis.  High  incident  angles  used  for  the  LT  AIN  sample  is  because  the  anisotropy  of  such  a  LT 
material  is  not  an  issue  in  the  material  quality  study. 

The  optical  constants  HT  and  LT  AIN  films  obtained  from  VASE  data  analysis  are  shown  in  Fig.  6.  The  HT  AIN  is  fully 
transparent  below  5  eV,  but  the  absorption  of  LT  AIN  starts  from  below  3  eV.  This  long  absorption  tail  below  the 
fundamental  band  gap  (-6.2  eV)  of  AIN  may  be  due  to  the  high-density  structure  defects  caused  by  both  low  temperature 
growth  and  large  lattice-mismatch  between  AIN  and  sapphire  substrate,  the  same  reason  caused  the  E2  phonon  broadening. 
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The  parametric  semiconductor  model  is  Kramer-Kronig  (K-K)  consistent,  the  K-K  fitting  is  shown  in  Fig.  7.  The  solid  lines 
are  the  best  £1  result  from  the  VASE  data  analysis,  and  the  dots  represent  the  calculated  Ei  from  the  known  e2  values  using  the 
K-K  relation.  Fig.  7  (a)  and  (b)  are  for  HT  and  LT  AIN,  respectively.  The  final  results  of  all  fitting  parameters  are  given  in 
Table  2.  The  K-K  relation  used  in  the  calculation  is17 


1  1  S  E?  -  (E)2 


9 

+  -P 
n  „ 


0.75eV 


e6-(eY 


For  a  specified  material,  the  K-K  integral  is  numerically  evaluated  to  calculate  £]  valufes  from  £2.  The  model  then  adds  the 


contribution  from  two  non-broadening  oscillators  and  a  fixed  offset  to  account  for  contribution  outside  the  experimental 
measuring  range.  E,  is  the  energy  location  of  a  non-broadening  oscillator  that  is  added  to  simulate  contribution  outside  the 
integration  range;  A,  is  the  magnitude  of  the  oscillator  located  at  E,  position.  In  our  calculation  one  oscillators  were  used 
outside  the  measured  region. 


Photon  Energy  (eV) 

Fig.  6.  Optical  constants  of  (a)  HT  AIN,  and  (b)  LT  AIN. 
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5.  Summary 


The  optical  anisotropy  nature  of  AIN  films  grown  on  c-plane  sapphire  substrates  was  revealed  by  Raman  scattering  and 
seneralized  variable  angle  spectroscopic  ellipsometry.  Optical  phonons  measured  from  a  1  |im  thick  HT  AIN  agrees  very 
well  with  the  theoretical  prediction.  Optical  constants  of  HT  and  LT  AIN  are  determined  by  TVASE  and  RVASE  along 
with  the  structure  parameters.  The  broadening  in  E?  phonon  and  long  absorption  tail  indicate  that  LT  AIN  has  very  high 
defect  density  due  to  low  temperature  growth  and  large  lattice-mismatch. 
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Standard  variable  angle  spectroscopic  ellipsometry  (VASE)  has  been  employed  to  study  the 
ordinary  optical  dielectric  response  of  hexagonal  gallium  nitride  (GaN)  thin  films — an  important 
material  for  blue  and  ultraviolet  light  emitting  device  applications.  The  GaN  films  were  grown  by 
molecular  beam  epitaxy  on  c-plane  sapphire  substrates  (a-Al203).  Room  temperature  isotropic  and 
anisotropic  mode  VASE  measurements  were  made  at  angles  of  incidence  between  of  20°  and  80°. 

Evidence  of  anisotropy  was  observed  from  the  anisotropic  mode  measurements,  reflecting  the  nature 
of  wurtzite  crystal  structure  of  GaN.  The  sizable  off-diagonal  elements  (Aps  and  Asp )  of  the  Jones 
matrix  indicate  that  the  optical  axis  ( c )  of  the  c-plane  sample  are  slightly  off  from  the  surface  normal 
due  to  a  small  miscut  of  substrates.  VASE  data  simulations  by  isotropic  and  anisotropic  models 
indicate  that  the  anisotropic  effect  on  both  diagonal  and  off-diagonal  elements  of  the  Jones  matrix 
can  be  minimized  to  a  negligible  level  at  small  angle  of  incidence.  Thus  the  ordinary  optical 
dielectric  functions  (£_L{c))  are  precisely  determined  by  the  isotropic  mode  VASE  measurements 
at  angles  of  incidence  between  20°  and  40°  in  the  range  of  0.75-6.5  eV.  The  VASE  data  were 
analyzed  by  a  model  dielectric  function  based  on  the  GaN  critical  point  structure,  which  allows  for 
a  nonzero  extinction  coefficient  k  below  the  band  gap.  The  thicknesses  of  these  GaN  films  are 
accurately  determined  via  the  analysis  as  well.  ©  2000  American  Institute  of  Physics. 
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I.  INTRODUCTION 

The  wide  band-gap  semiconductor  GaN  and  related  ma¬ 
terials,  with  their  excellent  thermal  conductivity,  large  break¬ 
down  field,  and  resistance  to  chemical  attack,  have  a  very 
promising  application  potential  for  both  high  temperature 
electronic  devices  and  short  wavelength  optical  emitters.1,2 
The  recent  development  of  high-brightness,  blue  and  green 
light  emitting  diodes  (LEDs),3  room  temperature  pulsed,4 
and  continuous-wave  quantum  well  lasers5  has  greatly  en¬ 
couraged  researchers  to  continue  the  work  on  these  materi¬ 
als.  It  is  essential  to  study  the  optical  properties  of  these 
optoelectronic  materials,  such  as  the  optical  dielectric  func¬ 
tions  (or  the  optical  constants)  of  GaN,  for  the  further  im¬ 
provement  of  device  performance.  The  optical  constants  of 
GaN  have  been  studied  intensively  since  the  mid-  to  late- 
1 960’s.  Typically,  the  ordinary  refractive  indices  (n)  were 
measured  in  the  visible  and  IR  spectral  range  (370-2000 
nm),  and  the  extinction  coefficients  ( k )  within  this  region 
were  either  ignored  or  assumed  as  zero,6-9  which  is  true  in 
theory  for  the  region  below  band  gap.  In  fact,  there  is  an 
absorption  tail  below  band  gap  for  most  GaN  films  due  to  the 
existence  of  defect  states  and  impurities.  Therefore,  it  would 
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be  more  accurate  if  a  nonzero  k  model  can  be  employed  for 
the  data  analysis.  The  recent  works  by  two  different  groups 
have  extended  the  index  of  refraction  data  to  the  region 
above  band  gap  using  spectroscopic  ellipsometiy.8,9  The 
nonzero  absorption  below  band  gap  was  also  observed  by 
other  groups.10,11 

Most  of  the  high  quality  GaN  films  were  grown  on 
c-plane  (0001)  sapphire  substrates  either  by  molecular  beam 
epitaxy  (MBE)  or  metalorganic  vapor  phase  epitaxy 
(MOVPE)  techniques.  Usually  GaN  grown  on  sapphire  sub¬ 
strates  has  a  wurtzite  structure  (a-GaN).  Nominally  the  op¬ 
tical  axis  (c)  of  the  GaN  film  is  perpendicular  to  the  sample 
surface,  and  parallel  to  the  (c)  of  sapphire.  The  fundamental 
band  gap  of  wurtzite  GaN  is  about  3.4  eV.7  As  the  growth 
techniques  are  being  improved,  higher  quality  GaN  films  are 
being  produced  with  lower  defect  densities,  lower  back¬ 
ground  electron  concentration,  and  higher  carrier  mobilities. 
More  precise  measurements  of  optical  properties  thus  are 
needed  for  future  research  and  fabrication  activities.  More¬ 
over,  GaN  is  an  anisotropic  crystal  (uniaxial),  and  this  an¬ 
isotropy  cannot  be  simply  neglected  during  the  optical  data 
acquisition  even  for  a  nominal  c-plane  sample.  The  existence 
of  optical  anisotropy  will  affect  both  the  diagonal  and  off- 
diagonal  elements  of  the  Jones  matrix.  If  the  optical  axis  is 
strictly  perpendicular  to  the  sample  surface  (perfect  c-plane 
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situation),  the  Jones  matrix  is  diagonal.  This  may  not  be  true 
for  most  GaN  films  grown  on  c-plane  sapphire  substrates, 
since  the  crystal  is  usually  cut  with  an  error  of  several  tenths 
of  degree  so  that  the  off-diagonal  elements  of  the  Jones  ma¬ 
trix  will  not  vanish  completely.  In  this  situation,  the  standard 
reflection  variable  spectroscopic  ellipsometry  (RVASE)  is 
not  reliable  because  a  portion  of  light  scattering  {p-s  or 
s-p)  has  been  neglected.  Here  “/?”  and  ‘V’  represent  the 
light  polarized  parallel  (p)  and  perpendicular  (5)  to  the  plane 
of  incidence,  respectively.  This  means  that  it  is  necessary  to 
use  the  generalized  ellipsometry  technique  to  detect  the  off- 
diagonal  elements  or  to  find  a  condition  at  which  the  stan¬ 
dard  ellipsometry  is  still  effective  with  negligible  off- 
diagonal  elements  Aps  and  Asp . 

It  would  be  ideal  and  straightforward  to  use  generalized 
ellipsometry  (GE)  to  determine  both  ordinary  and  extraordi¬ 
nary  dielectric  functions  of  GaN  simultaneously.  The  prob¬ 
lem  of  using  GE  on  a  slightly  misoriented  c-plane  GaN  is  the 
uncertainty  of  the  miscut  angle  of  the  sapphire  substrate  that 
makes  it  impossible  to  set  a  uniaxial  model  the  for  VASE 
data  analysis. 

Furthermore,  even  if  an  ideal  c-plane  GaN  is  available, 
the  standard  way  (such  as  near-Brewster  angle  of  incidence) 
to  perform  an  ellipsometric  measurement  for  an  isotropic 
sample  may  be  not  adequate  for  the  anisotropic  structure 
because  both  ellipsometric  parameters  'VT  and  A  (related  to 
the  diagonal  elements  of  Jones  matrix)  are  also  affected 
greatly  by  the  optical  anisotropy  (due  to  the  extraordinary 
response  for  p  polarization).  This  may  cause  a  few  percent 
error  in  the  resulted  optical  constants.  In  order  to  precisely 
measure  the  ordinary  optical  constants,  one  has  to  use  special 
incident  angles,  other  than  the  standard  approach  used  for 
isotropic  samples,  to  minimize  the  extraordinary  response 
(the  electric  field  Ell(c))  related  with  the  p  polarization  of 
the  incoming  light.  It  is  especially  important  to  obtain  highly 
accurate  ordinary'  dielectric  functions  for  the  determination 
of  the  extraordinary'  dielectric  functions  of  GaN,  since  the 
difference  between  them  is  just  a  few  percent  in  n  ratio.12 

The  ordinary  dielectric  functions  of  GaN  can  be  deter¬ 
mined  either  by  5-polarized  reflectance12  or  by  VASE  mea¬ 
surements.  In  theory  the  ^-polarized  reflectance  would  be 
ideal  for  the  determination  of  the  ordinary  dielectric  func¬ 
tions  of  a  uniaxial  anisotropic  sample,  since  the  polarization 
is  always  perpendicular  to  the  optical  axis  of  a  c-plane 
sample  regardless  of  the  angle  of  incidence.  The  ^-polarized 
reflection  coefficient  is  only  dependent  upon  the  ordinary 
dielectric  functions.12  But  in  practice,  ^-polarized  reflectance 
measurements  may  not  be  the  best  choice  since  they  do  not 
measure  the  phase  change.  For  most  of  the  samples,  the 
phase  information  is  necessary  and  critical  for  the  surface 
layer  analysis.  Without  the  phase  information,  the  optical 
properties  determined  by  the  reflectance  data  will  not  be  ac¬ 
curate.  It  is  especially  true  for  a  sample  with  some  native 
oxide  or  surface  roughness  on  top.  Another  major  drawback 
of  the  reflectance  technique  is  that  the  measured  data  are 
strongly  dependent  on  the  light  source  stability.  Any  tempo¬ 
ral  fluctuation  and  spatial  variation  may  cause  a  large  error  in 
the  reflectance  data.  On  the  other  hand,  ellipsometry  is  a  two 
parameter  ('R  and  A)  measurement  technique.  The  parameter 


A  contains  the  phase  information  needed  for  the  surface  layer 
uncertainty  determination.  Moreover  ellipsometric  data  will 
not  be  affected  by  the  light  source  instability.  One  problem 
for  using  VASE  to  determine  the  ordinary  dielectric  func¬ 
tions  is  the  extraordinary  response  caused  by  the  p-polarized 
light  at  nonzero  angles  of  incidence.  An  alternative  way  to 
reduce  this  extraordinary  response  is  to  use  a  small  angle  of 
incidence  in  the  VASE  measurement  to  keep  the  p-polarized 
part  nearly  perpendicular  to  the  optical  axis  (c). 

In  this  article,  the  optical  anisotropy  of  the  GaN  films 
grown  on  c-plane  sapphire  substrates  was  characterized  by 
generalized  variable  angle  spectroscopic  ellipsometry.  The 
impact  of  optical  anisotropy  on  the  determination  of  ordinary 
dielectric  functions  is  discussed  through  generated  reflectiv¬ 
ity  and  measured  VASE  data  by,  isotropic  and  anisotropic 
models.  Using  a  small  angle  of  incidence  in  the  VASE  data 
acquisition,  the  ordinary  dielectric  functions  of  GaN  are  pre¬ 
cisely  determined  via  a  parametric  semiconductor  model 
(with  nonzero  k)  in  a  wide  energy  range  between  0.75  and 
6.5  eV.  Furthermore,  consistency  with  Kramers -Kronig  con¬ 
ditions  is  maintained. 


II.  THEORY 

A.  Standard  ellipsometry 

Spectroscopic  ellipsometry  is  an  experimental  method 
that  can  be  used  to  measure  the  change  of  the  polarization 
state  after  light  is  reflected  from  or  transmitted  through  an 
optical  medium.  This  medium  can  be  a  bulk  substrate  or  a 
multiple  layer  structure.  By  analyzing  the  measured  ellip¬ 
sometry  parameters  if/  and  A  with  an  appropriate  physical 
model,  the  thickness  and  optical  dielectric  functions  of  each 
layer  can  be  determined.  The  parameters  ^and  A  are  related 
to  the  complex  ratio  of  reflection  (or  transmission)  coeffi¬ 
cients  for  light  polarized  parallel  (p)  and  perpendicular  (j)  to 
the  plane  of  incidence.13  For  isotropic  material  systems, 

P=^  =  tan((A)e'A.  (]) 

The  electric-field  reflection  coefficient  at  an  incident  angle  of 
(j>  is  defined  as  Rp  ( Rs )  for  p-  (j’)-polarized  light.  They  are 
the  diagonal  elements  of  the  Jones  matrix  in  a  nondepolar- 
ized  system: 


[*^]  sample 


0 


0 

Rs 


(2) 


The  parameters  if/  and  A  are  not  only  dependent  on  dielectric 
functions,  but  also  on  the  surface  condition,  sample  structure, 
and  other  properties  such  as  the  optical  anisotropy. 

The  pseudodielectric  function  (e)  can  be  obtained  di¬ 
rectly  from  the  measured  if/  and  A  values: 


<e)  =  <ei>-*"I(e2>  =  sin2  4> 


1  +  tan2  (j)  | 


(3) 


For  the  case  of  air  over  a  bare  bulk  material  with  a  perfectly 
smooth  surface,  the  pseudodielectric  and  intrinsic  dielectric 
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function  are  identical.  To  determine  the  optical  constants  of  a 
thin  film  on  a  substrate,  VASE  data  must  be  analyzed  using 
a  parametric  model  that  is  adjusted  to  fit  the  measured  data. 
A  regression  analysis  is  usually  used  to  vary  the  model  pa¬ 


rameters  (e.g.,  the  optical  constants  or  layer  thickness,  etc.) 
until  the  calculated  and  measured  values  match  as  closely  as 
possible.  This  is  done  by  minimizing  the  mean  square  error 
(MSE)  function,  defined  as: 


MSE= 


<P(h  Vi,<t>j)  ~  4f{h  Vii&j) 

2 

+ 

A(hVj, 4>j)  ~  &.c(hVj , <£/)! 2) 

»  N-MTj  1 

V i,jA  .  , 

(4) 


where  N  is  the  total  number  of  experimental  observations 
(for  ellipsometry  measurements,  there  are  two  observations 
and  A,  for  each  data  point),  M  is  the  number  of  fit  param¬ 
eters,  cr  is  the  measured  standard  deviation  of  the  measure¬ 
ment,  <$>  is  the  external  angle  of  incidence,  h  v  is  the  photon 
energy,  and  i  and  j  are  used  to  sum  over  all  the  photon 
energies  and  external  angles  of  incidence,  respectively. 

B.  Generalized  ellipsometry 

Generalized  ellipsometry  was  first  introduced  by  Azzam 
and  Bashara.14  The  recent  developments  make  this  technique 
more  capable  and  powerful.15,16  For  an  anisotropic  material 
system,  the  off-diagonal  elements  of  the  Jones  matrix  are  not 
necessary  to  be  zero.  In  the  reflection  (RVASE)  configura¬ 
tion, 


[^]  sample 


Rpp  R sp 
Rps  Rss 


(5) 


By  using  the  same  approach  with  considerably  more  compu¬ 
tation,  ij/and  A  values  can  still  be  determined,  but  the  Fourier 
coefficients  related  to  if/  and  A  become  more  complicated.17 
The  generalized  ellipsometric  parameters  are  defined  as  be¬ 
low: 


the  Fourier  coefficients14  are  the  functions  of  Aps  and  Asp , 
they  are  dependent  upon  the  incident  polarization  in  aniso¬ 
tropic  materials. 

III.  EXPERIMENT 

The  GaN  film  (about  1  p m  thick)  used  in  this  work  was 
grown  by  MBE  on  c-plane  sapphire  substrate.  The  cutting 
error  of  crystal  orientation  is  within  1°.  The  backside  of  the 
substrate  was  unpolished  which  is  sufficiently  rough  for  the 
reflection  VASE  and  polarized  reflectivity  measurements. 
The  RVASE  optical  measurements  of  both  isotropic  (stan¬ 
dard  ellipsometry)  and  anisotropic  modes  (generalized  ellip¬ 
sometry)  were  performed  in  the  energy  range  0.75-6.5  eV 
with  a  0.02  eV  increment  at  room  temperature.  In  order  to 
reveal  the  angular  dependence  of  the  optical  anisotropy  of  a 
c-plane  sample,  the  generalized  ellipsometry  measurements 
were  performed  at  1.4  eV  with  the  angle  of  incidence  rang¬ 
ing  from  20°  to  80°  with  1°  increments.  Silicon  and  InGaAs 
photodiodes  are  employed  in  the  VASE  system.  They  are 
polarization  state  insensitive,  and  linear  over  a  broad  range 
of  beam  intensity.  The  monochromator  was  placed  right  after 
the  light  source  for  wavelength  scanning. 

IV.  RESULTS  AND  DISCUSSION 


AnE=~j^~~ian  tJ/nEe‘AnE’  (6) 

Aps=TL^ tan  (Ap^'Api.  (7) 

nPP 

p 

Asp=-^L  =  tan  i Pspe‘*‘p,  (8) 

where  the  Aps  and  Asp  describe  how  much  amount  p-  or 
.y-polarized  light  becomes  s-  or  p-polarized  light  after  the 
reflection,  respectively.  The  AnE  is  the  diagonal  matrix  ele¬ 
ment. 

In  summary,  the  generalized  ellipsometry  is  a  technique, 
which  can  be  used  to  determine  all  the  elements  of  the  Jones 
matrix  for  arbitrarily  anisotropic  and  homogeneous  layered 
systems  with  nonscalar  dielectric  susceptibilities.  Unlike  the 
isotropic  mode  VASE  measurements,  the  anisotropic  mode 
measurements,  i.e.,  generalized  ellipsometry  are  made  with  a 
fixed  pattern  for  the  polarizer  position  (e.g.,  p  =  — 60°, 
-30°,  0°,  +30°,  +60°,  and  90°)  at  each  wavelength.  Since 


In  order  to  reveal  the  optical  anisotropy  of  GaN  on 
c-plane  sapphire,  a  series  of  anisotropic  mode  VASE  mea¬ 
surements  have  been  made.  Two  of  them  are  shown  in  Fig. 
1.  The  RVASE  data  in  Fig.  1(a)  were  taken  at  an  angle 
incidence  of  60°,  which  is  near  the  pseudo-Brewster  angle.  It 
can  be  seen  clearly  that  the  off-diagonal  elements  Aps  and 
Asp  are  not  negligible,  especially  for  Aps  below  band  gap. 
The  strong  Aps  fringes  appeared  below  3.4  eV  are  due  to  the 
thin  film  interference.  Notice  that  they  have  the  same  period 
as  the  AnE  does.  For  photon  energies  larger  than  band-gap 
value,  there  are  no  interference  fringes  since  the  GaN  film  is 
no  longer  transparent  in  this  region  (the  extinction  coefficient 
k  is  very  large).  Nonzero  Aps  and  Asp  indicate  that:  (1)  GaN 
film  has  strong  optical  anisotropy;  (2)  the  (c)  axis  of  GaN  is 
off  the  surface  normal.  For  such  a  sample,  the  Aps  and  Asp 
have  to  be  taken  into  account  in  the  data  analysis  unless 
special  measurement  conditions  are  employed  under  which 
the  anisotropic  effect  can  be  greatly  minimized.  The  nonzero 
Aps  and  Asp  means  that  there  is  cross-conversion  between  p- 
and  j-polarized  light.  We  cannot  set  up  an  anisotropic  model 
to  do  the  generalized-VASE  data  analysis  and  solve  both 
ordinary  and  extraordinary  dielectric  functions  simulta- 
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FIG.  1.  Anisotropic  mode  RVASE  data  of  c-plane  GaN/sapphire.  (a)  V  as  a 
function  of  photon  energy  at  60°  angle  of  incidence;  (b)  'vk  as  a  function  of 
incident  angle  at  1.4  eV. 


neously,  since  the  offset  angle  of  (c)  axis  from  surface  nor¬ 
mal  is  unknown.  An  easier  way  to  determine  the  dielectric 
functions  of  GaN  is  to  solve  them  separately,  and  to  con¬ 
struct  an  isotropic  model  instead  of  an  anisotropic  model  to 
solve  for  the  ordinary  dielectric  functions.  Therefore,  special 
measurement  conditions  must  be  used  to  suppress  the  Aps 
and  Asp  values,  and  to  allow  us  to  treat  the  anisotropic  film 
as  a  pseudoisotropic  material  without  losing  accuracy.  To 
find  such  measurement  conditions,  an  angular  dependent 
measurement  was  made  at  certain  wavelength  (1.4  eV),  as 
shown  in  Fig.  1(b).  The  angular  increment  is  1°.  The  off- 
diagonal  elements  are  detectable  in  the  angle  range  of  50°~ 
70°,  and  it  has  a  peak  value  at  about  60°,  which  is  near  the 
pseudo-Brewster  angle  at  1.4  eV.  We  have  measured  curves 
similar  to  that  in  Fig.  1(b)  at  other  photon  energies  with  the 
peak  positions  slightly  varied  due  to  different  Brewster 
angles  at  the  different  wavelengths.  In  the  photon  energy 
range  of  0.75-6.5  eV,  the  Aps  peak  position  can  change  from 
about  59°  to  67°,  i.e.,  Aps  may  be  detectable  from  approxi¬ 
mately  49°  to  77°  angle  of  incidence.  In  other  words,  the  Aps 
and  Asp  are  negligible  at  the  angle  of  incidence  below  40°  or 
above  80°  for  an  anisotropic  sample  with  slightly  (c)  axis 
offset.  The  above  measurement  conditions  are  only  useful  to 
minimize  the  off-diagonal  elements  ( Aps  and  Asp)  related  to 
the  misoriented  c-plane  GaN  film.  To  determine  ordinary 


Photon  Energy  (eV) 

FIG.  2.  Simulated  p-polarized  reflectivity  data  of  a  GaN  (2  /xm)/c-sapphire 
structure  by  isotropic  (solid  line)  and  anisotropic  (dash  line)  models  at  80° 
angle  of  incidence. 


dielectric  functions  of  GaN  more  precisely,  the  anisotropic 
effect  included  in  the  diagonal  elements  has  to  be  taken  into 
account  in  the  VASE  data  analysis. 

Usually  the  ordinary  dielectric  functions  of  anisotropic 
GaN  are  measured  by  standard  ellipsometry  at  large  angle  of 
incidence  such  as  60°  (Ref.  9)  and  70°, 18  and  modeled  as  an 
isotropic  material.  In  general,  results  from  these  studies  are 
suitable  as  provisional  data.  However,  they  lack  accuracy 
because  the  sample  misorientation  has  been  neglected  and 
because  the  anisotropic  nature  is  overlooked  in  the  data 
analysis.  This  may  become  a  problem  when  precise  determi¬ 
nation  of  the  extraordinary  dielectric  functions  is  required. 

To  illustrate  anisotropic  impact  on  the  diagonal  ele¬ 
ments,  the  difference  of  isotropic  and  anisotropic  optical  re¬ 
sponse  generated  from  ideal  c-plane  models  (no  misorienta¬ 
tion)  are  shown  in  Figs.  2  and  3.  Two  simulated 
GaN/c -sapphire  models  with  exactly  the  same  multilayered- 
structures  (2-yum  GaN  film  with  4%  thickness  nonunifor¬ 
mity)  were  constructed  for  the  data  simulation.  In  the  isotro¬ 
pic  model,  GaN  is  assumed  to  be  an  isotropic  material,  with 
ordinary  dielectric  functions  only.  In  the  anisotropic  model, 
GaN  is  assumed  to  be  an  anisotropic  material,  with  both 
ordinary  and  extraordinary  dielectric  functions.  The  ordinary 
and  extraordinary  functions  used  in  the  two  models  are  from 
a  study  on  GaN/m -sapphire. 19  The  average  difference  be¬ 
tween  extraordinary  and  ordinary  in  n  values  is  about  +4% 
in  the  range  of  0.75-6.5  eV.  Dielectric  functions  of  sapphire 
used  for  the  simulations  and  VASE  analysis  are  from  our 
previous  work.20  Data  generated  by  the  isotropic  and  aniso¬ 
tropic  models  are  identical  everywhere  because  the 
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FIG.  3.  Simulated  RVASE  data  of  GaN  (2  /xm)/c-sapphire  structure  by 
isotropic  (solid  line)  and  anisotropic  (dash  line)  models,  (a)  Large  angle  of 
incidence  (<£  =  80°);  (b)  small  angle  of  incidence  (<£  =  20°). 


5-polarized  light  can  only  cause  the  ordinary  dielectric  re¬ 
sponse  that  is  independent  of  the  angle  of  incidence.  But  for 
nonzero  degree  angle  of  incidence,  the  ^-polarized  incident 
light  can  cause  both  ordinary  and  extraordinary  responses  in 
an  anisotropic  material,  which  will  be  different  from  the  pure 
ordinary  response  in  an  isotropic  material,  as  shown  in  Fig. 
2.  The  reflectivity  difference  between  isotropic  and  aniso¬ 
tropic  model  is  about  3%  at  6  eV,  and  this  difference  is 
strongly  dependent  on  the  angle  of  incidence.  The  difference 
decreases  as  the  angle  of  incidence  decreases. 

Generated  VASE  data  with  different  angle  of  incidence 
for  the  two  models  are  shown  in  Fig.  3.  The  angle  of  inci¬ 
dence  is  80°  for  Fig.  3(a),  and  20°  for  Fig.  3(b).  Notice  that 
the  T'  difference  in  ratio  between  two  models 
([^aniso-^isoJ/^iso)  in  Fig-  3(a)  is  about  20  times  larger 
than  in  Fig.  3(b).  The  absolute  difference  between  values 
generated  by  two  models  in  Fig.  3(b)  is  close  to  the  measure¬ 
ment  error  bar  ±0.007,  which  is  within  the  sensitivity  limi¬ 
tation  of  the  instrument.  Therefore,  the  anisotropic  effect  on 
the  diagonal  elements  can  be  greatly  minimized  by  using  a 
small  angle  of  incidence.  In  addition,  the  small  angle  of  in¬ 
cidence  is  also  good  for  suppressing  the  s-p  or  p~s  cross 
conversion.  Comparing  with  the  measurement  error  bar,  we 
can  conclude  that  the  anisotropic  GaN  film  can  be  treated  as 
isotropic  material  at  small  angles  of  incidence  (0^40°). 

According  to  the  above  analysis,  VASE  data  at  a  small 
angle  of  incidence  (20°,  40°)  were  taken  from  a 
GaN/c-apphire  sample  for  the  determination  of  ordinary  di¬ 
electric  functions  of  GaN.  The  sample  structure  was  de¬ 
scribed  by  a  three-layer  model,  including  sapphire  substrate, 
an  isotropic  GaN  film  and  a  surface  layer  counting  for  the 


FIG.  4.  Structural  model  sketch  of  GaN/sapphire  with  a  surface  rough  layer. 


surface  roughness  of  GaN  film,  as  shown  in  Fig.  4.  The 
VASE  data  and  model  fittings  are  shown  in  Fig.  5.  The  top 
surface  layer  is  described  by  a  Bruggeman  effective  media 
approximation  (EMA)21,22  for  a  layer  consisting  of  GaN  and 
voids.  The  model  dielectric  function  (MDF)  of  GaN,  was 
represented  by  a  simplified  Adachi  model,18  which  consisted 
of  two  critical  points  (CP)  describing  the  interband  transi¬ 
tions  (E0  and  E}  gaps),  and  two  exciton  contributions  at  E0 
and  Ei  edges.  With  the  assumption  of  parabolic  band  shape 
for  both  valence  and  conduction  bands,  the  contribution  of 
E0  to  e(E)  is18,23 

e(0>(£)  =  A0E0-  15{*o2[2  -  ( 1  +  *o)°  5~  ( 1  -  ATo)0'5]}- 

(9) 
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FIG.  5.  RVASE  data  and  model  fittings  of  a  GaN/c -sapphire  sample,  (a)  V 
data  and  the  best  model  fit  at  20°  and  40°  angles  of  incidence;  (b)  A  data  and 
the  best  model  fit  at  20°  and  40°  angles  of  incidence.  The  solid  line  repre¬ 
sents  the  best  model  fit.  The  circle  and  square  represent  the  data  taken  at  20° 
and  40°  angles  of  incidence,  respectively. 
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FIG.  6.  Ordinary  dielectric  functions  of  GaN  obtained  from  the  small  angle 
of  incidence  RVASE  data  analysis  of  a  GaN/ c -sapphire  sample. 
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FIG.  7.  Kramers-Kronig  consistency  check  of  the  ordinary  dielectric  func¬ 
tions  of  GaN. 


The  Ex  contribution  to  e(£)  can  be  expressed  as 

e<15(£,)=i41^[’2  ln(l  —  Xi)>  (10) 

and 

Xj  =  {E+iTj)IEj,  (11) 

where  j  can  be  0  or  1  for  Eqs.  (9)  and  (10),  respectively.  Aj , 
Ej ,  and  F;  are  amplitude,  transition  energy,  and  broadening 
parameter  of  the  E0  and  Ex  CP  structures,  respectively.  The 
contributions  to  e(£)  from  excitons  at  £0  and  Ex  edges  are 
approximated  by  two  discrete  damped  Lorentzian  line  shapes 
with  complex  amplitudes: 

Ajx^MiPjx) 

E  -E~i T  ’ 

where  j  =  0,  1  for  £0  and  Ex ,  respectively.  Ejx ,  T jx ,  and  Ajx 
are  center  energy,  broadening,  and  amplitude  of  each  exci- 
ton,  respectively,  and  p]X  is  the  phase  factor  of  the  complex 
amplitude.  The  total  dielectric  response  due  to  interband 
transitions  and  excitonic  contributions,  in  the  energy  range  of 
0.75-6.5  eV,  is  then  given  by 

e(E)  =  ex+  e(0)(£)  +  e(0j:)(£)  +  e(1)(£)+  e(1*>(£). 

(13) 

The  dielectric  constant  accommodates  contributions  from 
CP  structures  at  energies  higher  than  Ex ,  which  are  not  di¬ 
rectly  included  in  the  model.  The  upper  CPs  are  not  consid¬ 
ered  in  this  model  is  due  to  the  limit  of  experimental  data  at 
6.5  eV.  This  parametric  dispersion  model  can,  without  com¬ 
promise  to  the  quality  of  data  fitting,  replace  the  usual  tabu¬ 
lated  optical  constants  lists  by  using  a  reasonable  small  set  of 
adjustable  parameters  (15  parameters  total  in  this  model  ex¬ 


cept  film  thickness  and  surface  parameters).  It  reduces  pa¬ 
rameter  correlation  in  data  fittings,  and  enforces  continuity 
on  the  optical  constants. 

With  a  small  amount  of  surface  roughness  (about  3.7  nm 
thick  EMA  layer  with  41%  voids  resulted  from  the  best  fit), 
the  VASE  data  at  both  angles  of  incidence  were  fitted  ex¬ 
tremely  well.  The  GaN  film  thickness  resulting  from  VASE 
data  analysis  is  about  1024  nm,  which  is  close  to  the  nominal 
value  of  1  /xm. 

The  ordinary  dielectric  functions  of  GaN  film  grown  on 
c-plane  sapphire  substrate  determined  in  this  work  are  shown 
in  Fig.  6.  The  fitting  parameters  of  GaN  ordinary  dielectric 
functions  are  listed  in  Table  I.  The  starting  values  of  the 
fitting  parameters  are  set  based  on  the  critical  point  structure 
of  GaN.18  Notice  that  the  £0  and  Ex  values  are  slightly  dif¬ 
ferent  from  those  given  in  Ref.  18.  This  may  be  due  to  the 
simplification  of  the  present  model  and  the  limited  spectral 
range  of  the  experimental  data.  The  sharp  turn  of  e2  near  the 
band  edge  and  the  very  weak  absorption  tail  indicate  that  this 
film  has  high  crystal  quality.  In  order  to  check  the  Kramer- 
Kronig  (KK)  consistency  of  the  resulting  GaN  ordinary  di¬ 
electric  functions,  a  KK  fitting  is  shown  in  Fig.  7.  The  circles 
are  the  best  ex  result,  and  the  solid  line  is  the  calculated  ex 
from  the  known  e2  values  using  the  KK  relation.  The  final 
results  of  all  fitting  parameters  are  given  in  Table  II.  The  KK 
relation  used  in  the  calculation  is24,25 


£f(£)  =  6fset+ 


Y  Ai 

&  E2-(E)2 


2  f6.5ev  £'e”eas(£') 

+  -p  >  dE' 

T*  J  0.75  eV  £  \tL) 


(14) 


TABLE  I.  Best  fitting  parameters  for  GaN  ordinary  dielectric  functions. 


CP  contribution  Aj(AJX) 

E,(E,X) 

Tj(rjx) 

Pjx 

£<°> 

30.15  (eV)15 

3.44  (eV) 

0.22  (eV) 

£«W 

1.07  (eV) 

3.42  (eV) 

0.19  (eV) 

-1.52 

2.90  6.94  (eV)  0.047  (eV) 

3.34  (eV)  6.92  (eV)  1.86  (eV) 
0.23 


TABLE  II.  Kramers-Kronig  fitting  parameters  for  GaN  ordinary  dielectric 
functions. 

Fitting  parameters 

Pole  No.  Et  A i  £jffsel 


-0.025 


1 

2 


7.491 

14.092 


2.001 

24.6 


1.04 
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For  a  specified  material,  the  KK  integral  is  numerically 
evaluated  to  calculate  ex  values  from  e2*  The  model  then 
adds  the  contribution  from  two  nonbroadening  oscillators 
and  a  fixed  offset  to  account  for  absorption  outside  the  ex¬ 
perimental  measurement  range.  E{  is  the  energy  location  of  a 
nonbroadening  oscillator  that  is  added  to  simulate  absorption 
outside  the  integration  range,  At-  is  the  magnitude  of  the  os¬ 
cillator  located  at  Ex  position.  In  our  calculation  two  oscilla¬ 
tors  were  used  outside  the  measured  region. 

V.  SUMMARY 

The  optical  anisotropy  of  GaN  films  grown  on  c-plane 
sapphire  substrates  was  revealed  by  the  generalized  variable 
angle  spectroscopic  ellipsometry.  This  anisotropy  can  be  ob¬ 
served  via  the  nonzero  off-diagonal  Jones  matrix.  These  off- 
diagonal  elements  are  nonnegligible  in  a  range  of  incident 
angles  for  slightly  miscut  samples.  However,  these  elements 
may  become  undetectable  beyond  this  range  of  angles.  At  a 
small  angle  of  incidence  (<£<40),  the  optical  anisotropic 
effect  on  both  off-diagonal  and  diagonal  elements  of  the 
Jones  matrix  can  be  greatly  reduced.  Thus  the  ordinary  di¬ 
electric  functions  of  GaN  are  accurately  determined  in  the 
range  of  0.75-6.5  eV  by  fitting  the  RVASE  data  with  an 
isotropic  CP  model  dielectric  function.  The  comparison  of 
the  calculated  and  the  experimental  ex  indicates  that  the 
dielectric  functions  obtained  in  this  study  are  KK  consistent. 

ACKNOWLEDGMENT 

This  work  was  supported  by  U.S.  Army  Research  Office 
under  Contract  No.  DAAG55-98-1-0462. 

1  H.  Morko?,  S.  Strite,  G.  B.  Gao,  M.  E.  Lin,  B.  Sverdlov,  and  M.  Bums,  J. 

Appl.  Phys.  76,  1363  (1994). 


2S.  Nakamura,  T.  Mukai,  and  M.  Senoh,  Appl.  Phys.  Lett.  64,  1687  (1994). 
3S.  Nakamura,  M.  Senoh,  N.  Iwasa,  and  S.  Nagahama,  Jpn.  J.  Appl.  Phys., 
Part  2  34,  L797  (1995). 

4S.  Nakamura,  M.  Senoh,  S.  Nagahama,  N.  Iwasa,  T.  Yamada,  T.  Mat¬ 
sushita,  H.  Kiyoku,  and  Y.  Sugimoto,  Jpn.  J.  Appl.  Phys.,  Part  2  35,  L74 
(1996). 

5S.  Nakamura,  M.  Senoh,  S.  Nagahama,  N.  Iwasa,  T.  Yamada,  T.  Mat¬ 
sushita,  H.  Kiyoku,  Y.  Sugimoto,  T.  Kozaki,  H.  Umemoto,  M.  Sano.  and 
K.  Chocho,  Appl.  Phys.  Lett.  72,  2014  (1998). 

6E.  Ejder,  Phys.  Status  Solidi  A  6,  445  (1971). 

7S.  Logothetidis,  J.  Petalas,  M.  Cardona,  and  T.  D.  Moustakas,  Phys.  Rev. 
B  50,  18017  (1994). 

8H.  Amano,  N.  Watanabe,  N.  Koide,  and  I.  Akasaki,  Jpn.  J.  Appl.  Phys., 
Part  2  32,  L1000  (1993). 

9G.  Yu,  G.  Wang,  H.  Ishikawa,  M.  Umeno,  T.  Soga,  T.  Egawa,  J.  Wa¬ 
tanabe,  and  T.  Jimbo,  Appl.  Phys.  Lett.  70,  3209  (1997). 

10W.  R.  L.  Lambrecht,  B.  Segall,  J.  Rife,  W.  R.  Hunter,  and  D.  K.  Wick- 
enden,  Phys.  Rev.  B  51,  13516  (1995).  - 
11  J.  Petalas,  S.  Logothetidis,  S.  Boultadakis,  M.  Alouani,  and  J.  M.  Wills, 
Phys.  Rev.  B  52,  8082  (1995). 

12G.  Yu,  H.  Ishikawa,  T.  Egawa,  T.  Soga,  J.  Watanabe,  T.  Jimbo,  and  M. 

Umeno,  Jpn.  J.  Appl.  Phys.,  Part  2  36,  LI 029  (1997). 

13 R.  M.  A.  Azzam  and  N.  M.  Bashara,  Ellipsometry  and  Polarized  Light 
(North-Holland,  Amsterdam,  1977). 

14 R.  M.  A.  Azzam  and  N.  M.  Bashara,  J.  Opt.  Soc.  Am.  62,  1521  (1972). 
15D.  W.  Berreman,  J.  Opt.  Soc.  Am.  62,  502  (1972). 

16 M.  Schubert,  Phys.  Rev.  B  53,  4265  (1996). 

17 M.  Schubert,  B.  Rheinlander,  J.  A.  Wollam,  B.  Johs,  and  C.  M.  Herzinger, 
J.  Opt.  Soc.  Am.  A  13,  875  (1996). 

18 T.  Kawashima,  H.  Yoshikawa,  S.  Adachi,  S.  Fuke,  and  K.  Ohtsuka,  J. 

Appl.  Phys.  82,  3528  (1997). 

19C.  H.  Yan  and  H.  Yao  (unpublished). 

20  H.  Yao  and  C.  H.  Yan,  J.  Appl.  Phys.  85,  6717  (1999). 

21 D.  E.  Aspnes,  J.  B.  Theeten,  and  F.  Hottier,  Phys.  Rev.  B  20,  3292  (1979). 
22H.  Yao,  B.  Johs,  and  R.  B.  James,  Phys.  Rev.  B  56,  9414  (1997). 

23  M.  Schubert,  J.  A.  Woollam,  G.  Leibiger,  B.  Rheinlander,  I.  Pietzonka,  T. 

Sass,  and  V.  Gottschalch,  J.  Appl.  Phys.  86,  2025  (1999). 

24H.  Yao,  B.  Johs,  and  R.  B.  James,  Phys.  Rev.  B  56,  9414  (1997). 

25 H.  W.  Yao,  J.  C.  Erickson,  H.  B.  Barber,  and  R.  B.  James,  J.  Electron. 
Mater.  28,  760  (1999). 


Manuscript  for  submitting  to  Journal  of  Applied  Physics 


Extraordinary  optical  dielectric  functions  of  anisotropic  hexagonal  GaN  film 
determined  by  variable  angle  spectroscopic  ellipsometry 

C.  H.  Yan,  H.  Yaoa), 

Center  for  Microelectronic  and  Optical  Material  Research,  and  Department  of  Electrical 
Engineering,  University  of  Nebraska,  Lincoln,  NE  68588-0511 

J.  M.  Yan  Hove,  A.  M.  Wowchak,  and  P.  P.  Chow 
SVT  Associates,  Inc 

J.  M.  Zavada 

US  Army  European  Research  Office,  London,  UK 

Abstract 

Variable  angle  spectroscopic  ellipsometry  (VASE)  has  been  employed  to  study 
the  extraordinary  optical  dielectric  response  of  hexagonal  gallium  nitride  (GaN)  thin  film 
grown  on  c-plane  sapphire  substrates  (a-Al203)  by  molecular  beam  epitaxy  (MBE). 

Room  temperature  VASE  measurements  were  made,  in  the  range  of  0.75  to  6.5  eV,  at  the 
angle  of  incidence  in  between  of  60  and  80  degree.  VASE  data  simulations  by  isotropic 
and  anisotropic  models  indicate  that  the  anisotropic  effect  can  be  detected  most- 
sensitively  at  large  angles  of  incidence  (near  the  pseudo-Brewster  angle).  Thus  the 
extraordinary  optical  dielectric  functions  (E  ||  <c>)  are  precisely  determined  by  standard 
VASE  measurements  at  60°,  70°,  and  80°  angles  of  incidence  in  the  range  of  0.75  to  6.5 
eV.  The  surface  roughness  was  also  considered  in  the  VASE  data  analysis  in  order  to 
separate  it  from  the  anisotropic  effect.  The  VASE  data  is  analyzed  by  a  model  dielectric 
function  based  on  the  GaN  critical  point  structure,  which  allows  for  a  non-zero  extinction 
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coefficient  k  below  the  band  gap.  The  thicknesses  of  these  GaN  films  are  accurately 
determined  via  the  analysis  as  well. 
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1.  Introduction 


The  wide  band  gap  semiconductor  GaN  is  an  important  material  for  light  emitting 
device  applications  in  the  green,  blue  and  UV  regions1, 2’ 3.  GaN  films  grown  on  sapphire 
substrates  usually  have  hexagonal  structure,  which  is  optically  anisotropic.  The  ordinary 
optical  dielectric  functions  of  GaN  have  been  studied  intensively  since  1960’s4, 5’ 6.  The 
most  recent  data  of  ordinary  dielectric  functions  were  determined  by  variable  angle 
spectroscopic  ellipsometry  (VASE)  with7  or  without8  the  consideration  of  surface 
roughness  and  the  optical  anisotropy.  But  the  extraordinary  dielectric  functions  of  GaN, 
on  the  other  hand,  are  not  well  documented.  Yu9  reported  the  extraordinary  functions 
below  the  fundamental  band  gap  measured  by  polarized  reflectance  without  the 
consideration  of  the  surface  roughness.  Polarized  reflectance  spectroscopy  is  a 
straightforward  technique  to  measure  the  extraordinary  response  of  GaN.  However  it 
may  be  lack  of  accuracy  because  not  only  reflectance  is  an  intensity  sensitive  method,  but 
also  it  is  lack  of  sensitivity  on  surface  roughness.  The  later  is  because  of  the  phase 
information  cannot  be  obtained  in  a  reflectance  spectrum.  But  this  does  not  mean  that  the 
reflectance  spectrum  will  not  be  affected  by  the  surface  roughness. 

In  principle,  the  best  way  to  determine  the  extraordinary  dielectric  functions  is 
using  the  generalized  ellipsometry  (GE),  which  can  accurately  distinguish  the  difference 
between  the  ordinary  and  extraordinary  dielectric  response  by  measuring  the  off-diagonal 
elements  of  Jones  matrix  (Aps  and  Asp)7.  But  the  GE  measurements  for  the  determination 
of  the  extraordinary  functions  require  non  c-plane  samples,  or  say  the  optical  axis  of  a 
GaN  film  has  to  be  certain  degree  away  from  the  surface  normal,  and  this  angle  has  to  be 
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well  determined  before  the  GE  measurements.  In  fact,  Most  of  high  quality  GaN  films 
were  grown  on  c-plane  (0001)  sapphire  substrates  either  by  molecular  beam  epitaxy 
(MBE)  or  metalorganic  vapor  phase  epitaxy  (MOVPE)  techniques.  The  optical  axis  <c> 
of  the  GaN  film  grown  on  c-plane  sapphire  is  perpendicular  to  the  sample  surface 
nominally,  and  parallel  to  the  <c>  of  sapphire.  In  this  case  the  off-diagonal  elements  of 
Jones  matrix  are  zero.  Therefore,  Measuring  a  GaN  film  grown  on  c-plane  sapphire  by 
GE  will  not  bring  any  information  about  the  off-diagonal  elements.  An  alternative  way  to 
determine  the  extraordinary  dielectric  functions  of  GaN  films  grown  on  c-plane  sapphire 
substrates  is  to  employ  the  standard  ellipsometry  under  some  special  conditions. 

For  an  anisotropic  material,  both  'F  and  A  will  be  affected  greatly  by  the  optical 

anisotropy.  This  effect  is  also  incident-angle  dependent.  At  certain  angle  of  incidence 
(near  pseudo-Brewster  angle),  extraordinary  optical  response  (E  ||  <c>)  has  the  maximum 
sensitivity.  At  small  angle  of  incidence,  the  extraordinary  response  has  almost  no 
sensitivity  in  the  VF  and  A.  This  is  how  the  ordinary  dielectric  functions  were  determined 

in  our  previous  study7.  In  this  paper,  the  extraordinary  dielectric  functions  of  a  GaN  film 
grown  on  c-plane  sapphire  substrates  were  determined  by  standard  variable  angle 
spectroscopic  ellipsometry  at  large  angles  of  incidence,  via  a  parametric  semiconductor 
model  (allow  nonzero-k)  in  a  wide  energy  range  of  0.75eV  to  6.5  eV,  and  they  are 
Kramers-Kronig  consistent.  The  impact  of  surface  roughness  on  the  VASE  data  analysis 
is  also  discussed  through  model-generated  reflectivity  and  VASE  data. 


4 


2.  Theory 


2.1.  Ellipsometry  for  isotropic  films 

For  an  isotropic  material  system,  the  ellipsometric  parameters  'F  and  A  are 
defined  as  below. 

p  =  tan(v^)e'A  (1) 

Ks 


Rs  and  Rp  are  the  total  reflection  coefficients  of  a  multiple-layer  structure  for  s  and  p 
polarized  light,  respectively.  They  are  the  diagonal  elements  of  Jones  matrix  in  a  no 
depolarized  system.  In  the  case  of  a  substrate  with  one  surface  overlayer,  light  can  see 
two  interfaces  (air-film,  and  film-substrate)  for  both  reflection  and  transmission.  The 
total  reflection  coefficients  (Rs  and  Rp)  are  expressed  as10: 


R  =r?2+r{3exp(-j2P) 
1  +  ri2r23  exP (rfiP) 


(2) 


R, 


=  rn  +  r23  exP 

1+r12r23  exP (~JW) 


(3) 


Where 


P  =  j_)N2  COS02 


(4) 


rn  5  r23 ,  rn  ,  and  r23  are  the  Fresnel  reflection  coefficients  of  a  single  interface.  The 


general  formats  of  p-  and  s-polarized  reflection  coefficients  are 
_  N J  COS0,  -N2COS<j>2 


ns2 


yv,  cos<p|  +yV2cos</)2 


(5) 
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(6) 


P  _  N2  COS0j  -  Ny  COS02 
12  Af2  COS0!  +  iVj  COS02 

N,  and  N2  are  the  complex  indices  of  reflection  of  the  upper  and  lower  medium.  P\  is 

the  angle  of  incidence,  and  02  is  the  refraction  angle  (they  are  related  by  Snell’s  law). 

The  V}/  and  A  are  not  only  dependent  on  dielectric  functions,  but  also  on  the  angle 

of  incidence,  surface  condition,  sample  structure,  and  other  properties  such  as  the  optical 

anisotropy  for  anisotropic  materials. 


2.2.  Ellipsometry  of  anisotropic  films 

For  a  uniaxial  anisotropic  film,  when  the  optical  axis  <0  is  perpendicular  to  the 
film  surface,  i.e.,  the  c-plane  case,  the  three-phase  total  reflection  coefficients  for  s-  and  p- 
polarized  light  are11 


R  =  rl2+r23exP  Hips) 
1 +  r\2r23  exp(-y'2  Ps) 

R  =r{2+r&exp(-j2pp) 
1  +  r12  r23  exp(-/2  Pp  ) 


Where  ri2  ,  r23  are  the  Fresnel  reflection  coefficients  at  the  first  and  second  interfaces  for 

P  P 

s-polarization.  rn  and  r23  are  the  Fresnel  reflection  coefficients  at  the  first  and  second 

interfaces  for  p-polarization.  They  are 
s  _Ny  cos0j  -N2±  cos02 

P\2  ~  _ _ c* _  (9) 

N\  COS01  +  N2i_  COS02 


r23  " 


N21  COS0]  -  COS 03 
N2\  COS 0i  +  A^3j_  COS02 


(10) 
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(11) 


rp  _  NqPis.  cos0i  -^i[(^2||)2  -(^i)2sin2^]^“ 
12  +  ^i[(^2||)2  -(^i)2sin2^]^“ 

p  _  ^3[(^2||)2  - A^i  sin2  </>,  ]  ^"-JV 2j|iv21cos03 
23  jV3[(^2||)2  ~  ^isin2  <Px]r+N^N2i_  cos(j)3 


Ps  and  Pp  represent  the  phase  difference  for  the  s-polarized  wave  and  p-polarized  wave. 


_2^[(Ar21)2-(Af1)2sin2«^1]- 
Ps  - j - 

_2^21[(^|)2-(iV1)2sin2^1]- 


(13) 

(14) 


Where  A  is  the  wavelength  of  the  monochromatic  light  in  vacuum  and  d  is  the  thickness  of 
the  anisotropic  film. 

Look  into  the  definition  of  ¥  and  A,  it  is  obvious  that  these  two  parameters 
depend  on  nzl  and  ^2|| ,  and  the  angle  of  incidence. 


3.  Experiments 

The  GaN  film  (about  lpm  thick)  studied  in  this  work  was  grown  by  MBE  on  c- 

plane  sapphire  substrate.  The  cutting  error  of  crystal  orientation  is  within  about  1°.  The 
backside  of  the  substrate  was  unpolished  which  is  rough  enough  for  making  good 
reflection  VASE  and  polarized  reflectivity  measurements  from  the  front  surface.  The 
RVASE  optical  measurements  of  both  isotropic  and  anisotropic  modes  were  performed  in 
the  energy  range  of  0.75  eV  to  6.5  eV  with  a  0.02  eV  increment  at  room  temperature. 
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Multiple  angles  of  incidence  ranging  from  60°  to  80°  were  employed  in  order  to  resolve 
the  optical  anisotropy  of  GaN  film  and  its  surface  roughness. 

4.  Results  and  discussion 

In  order  to  explore  the  extraordinary  response  of  GaN  on  c-plane  sapphire,  a 
optimum  measurement  condition  has  to  be  determined  before  the  data  taking.  Generated 
'F  values  at  an  arbitrary  photon  energy  (2  eV)  as  a  function  of  incident  angle  are  shown  in 

Fig  1 .  Data  in  Fig.  1  (a)  is  obtained  from  a  simulated  sample  structure  consisted  of  a  2  |im 
sapphire  film  and  a  silicon  substrate.  Fig.  1  (b)  is  a  similar  simulation  from  a  GaN  (2 
|im)/c-sapphire  structure.  Isotropic  and  anisotropic  model  mean  that  the  sapphire  and 

GaN  films  are  assumed  isotropic  or  anisotropic  in  data  generation.  The  ordinary  and 
extraordinary  dielectric  functions  of  sapphire  and  GaN  used  in  the  calculation  are  from 
our  previous  study7,  n.  It  is  obvious  that  'F  values  are  strongly  dependent  on  the  angle  of 

incidence  for  both  isotropic  and  anisotropic  models.  The  difference  between  isotropic 
and  anisotropic  model  is  significant  at  large  angles  near  the  pseudo-Brewster  angle  in  both 
cases.  At  low  incident  angles  (<50°),  the  difference  is  negligible.  Therefore,  60°,  70°  and 
80°  are  the  good  angles  for  resolving  the  extraordinary  response  of  the  GaN  film. 

For  a  real  sample,  surface  over-layer  is  always  an  important  issue  in  an 
ellipsometric  study,  since  the  native  oxide  or/and  surface  roughness  are  existed  on  many 
materials.  Surface  over-layer  can  alter  both  'F  and  A  significantly  and  will  affect  the 
values  of  the  extracted  dielectric  functions  from  the  VASE  data.  In  order  to  see  the  over- 
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layer  effect  on  a  GaN/sapphire  sample,  a  simulated  sample  is  constructed  which  consists 
of  an  EMA  (effective  media  approximation)  layer13  (50%  voids,  and  different  thickness), 
an  isotropic  GaN  film  (1  (im  thick),  and  a  c-plane  sapphire  substrate,  as  shown  in  Fig.  2. 
The  generated  VASE  data  as  a  function  of  incident  angle  is  shown  in  Fig.  3.  The  photon 
energy  in  Fig.  3(a)  and  (b)  is  5  eV  (above  the  GaN  band  gap),  and  2  eV  (below  the  band 
gap)  in  Fig.  3  (c)  and  (d).  It  is  obvious  that  both  'F  and  A  are  very  sensitive  to  the 
thickness  of  surface  over-layer,  and  the  most  sensitive  region  happens  near  the  Brewster 
angle,  for  both  5  eV  and  2  eV  cases.  Notice  that  A  is  much  more  sensitive  to  surface- 

roughness  than  VF,  and  the  sensitivity  could  be  in  one  monolayer  scale.  Therefore  the 

surface  roughness  can  be  accurately  determined  by  VASE.  This  is  one  of  the  main 
reasons  that  VASE  is  a  more  sensitive  measurement  than  reflectance  spectroscopy.  The 
surface  effect  can  also  be  observed  in  the  full  range  'F  and  A  spectra,  as  shown  in  Fig.  4. 
With  only  4  nm  difference  in  over-layer  thickness,  A  alters  significantly  at  both  angles  in 
the  above  band  gap  region.  It  is  also  shown  that  the  surface  effect  for  the  below  band  gap 
region  is  much  weaker. 

Generated  reflectance  data  with  and  without  surface  roughness  is  shown  in  Fig.  5. 
As  mentioned  above,  the  surface  effect  is  more  sensitive  in  the  above  band  gap  region  in 
the  reflectance  spectra.  It  is  also  sensitive  to  rough  surface  at  large  incident  angles 
(65°«}><850  in  this  GaN/sapphire  case).  Therefore  using  75°  as  the  angle  of  incidence9 
without  considering  the  surface  roughness  may  result  a  large  error  bar  in  the  resulted 
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dielectric  functions.  In  stead,  60°  may  be  a  better  choice  as  the  angle  of  incidence  if  using 
polarized  reflectance  spectroscopy  for  GaN  study. 

Based  on  the  above  discussions  we  know  that  both  anisotropy  and  surface 
roughness  affect  'P  and  A  spectra.  In  general,  surface  roughness  affects  A  more  than  'F, 
especially  at  the  near  Brewster  angle  region  and  above  band  gap  region.  While  optical 
anisotropy  affects  'F  more  than  A,  and  this  happens  in  the  whole  spectral  range,  and  high 

angles  of  incidence.  There  are  some  conditions  at  which  these  two  effects  may  be 
separated,  such  as  at  60°  angle  of  incidence  for  the  GaN/c-sapphire  case.  It  can  be  found 
in  Fig.  4  and  Fig.  5  that  there  is  a  very  weak  surface  roughness  effect  at  60°  angle  of 
incidence  while  the  anisotropy  effect  is  still  manageable  as  shown  in  Fig.  1  (b).  In  order 
to  resolve  both  the  extraordinary  response  and  the  surface  roughness,  a  series  of  VASE 
measurements  were  performed  at  large  angles  of  incidence,  i.e.,  60°,  70°,  and  80°.  The 
VASE  data  is  analyzed  with  a  three-phase  model  as  shown  in  Fig.  6.  The  fitting 
parameters  are  the  surface  over-layer  thickness  di,  its  void  percentage,  GaN  film 
thickness  d2  and  the  MDF  parameters  for  GaN  dielectric  functions14,15  of  the 
extraordinary  part.  The  GaN  ordinary  dielectric  functions7  and  sapphire  dielectric 
functions12  are  from  our  previous  study.  Some  of  the  final  fitting  results  are  shown  in  the 
Fig.  6.  The  GaN  film  is  1024  nm  thick,  which  is  closed  to  the  nominal  value  (1  pm).  The 

surface  over-layer  is  about  3.7  nm  thick  with  41%  voids.  The  VASE  data  and  the  best  fit 
are  shown  in  Fig.  7.  The  fitting  parameters  of  GaN  extraordinary  dielectric  functions  are 
listed  in  Table  4.1.  The  details  of  this  parametric  model  description  and  the  method  of 
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use  are  same  as  the  ordinary  part.  The  resulting  ordinary  and  extraordinary  dielectric 


functions  are  plotted  in  Fig.  8.  is  larger  than  ei±  everywhere  in  the  whole  measured 

range,  and  the  average  difference  between  them  is  about  3.5%.  The  abrupt  step  around 
band  gap  in  £ 2  curves  indicates  the  high  quality  of  this  film,  which  indicates  low  defect 
density  and  low  impurity  concentration. 

Kramer-Kronig  consistency  check  was  performed  on  both  ordinary  and 
extraordinary  dielectric  functions.  The  k-k  results  are  shown  in  Fig.  9,  and  Table  2  and  3. 
The  solid  dots  circles  are  the  best  6]  results  for  ordinary  and  extraordinary,  respectively. 


The  solid  line  is  the  calculated  Ei  from  the  known  e2  values  using  the  K-K  relation.  The 
K-K  relation  used  in  the  calculation  is16, 17 


e?(E)  =  £,o/m  + 


Ai 


.+ 


6-5eK  771  „meas  /  771  \ 

Ee;  ,(£V 

£'z-(£)2 


(15) 


0.75eV 


For  a  set  of  known  dielectric  functions,  £1  can  be  calculated  from  e2  using  the  above  K-K 
integral  in  the  measured  energy  region  (0.75  eV  to  6.5  eV).  The  model  then  adds  the 
contribution  from  two  non-broadening  oscillators  and  a  fixed  offset  to  account  for 
absorption  outside  the  experimental  measuring  range.  Et  is  the  energy  location  of  a  non¬ 
broadening  oscillator  that  is  added  to  simulate  the  absorption  outside  the  integration 
range.  At  is  the  magnitude  of  the  oscillator  located  at  E,  position.  In  our  calculation  two 
oscillators  were  used  outside  the  measured  region. 
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Summary 


Extraordinary  dielectric  functions  of  a  GaN  film  grown  on  c-plane  sapphire 
substrates  were  revealed  by  the  standard  variable  angle  spectroscopic  ellipsometry  in  the 
range  of  0.75  to  6.5  eV.  Model  simulations  indicate  that  large  incident  angles  (near  the 
pseudo-Brewster  angle)  have  the  best  sensitivity  on  resolving  the  extraordinary  response 
and  the  surface  roughness  simultaneously.  It  also  shows  that  VASE  is  a  better  choice 
than  the  reflectance  spectroscopy  when  surface  over-layer  exists.  The  real  part  of  the 
extraordinary  function  is  larger  than  that  of  the  ordinary  by  an  average  value  of  3.5%. 
The  comparison  of  the  calculated  £i  and  the  experimental  £j  indicates  that  the  dielectric 
functions  obtained  in  this  study  are  K-K  consistent. 
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Figure  captions 

Fig.  1 .  Generated  VASE  data  as  a  function  of  angle  of  incidence  at  2  eV.  (a)  Sapphire  (2 
pm  thick)  on  silicon  substrate,  (b)  GaN  (2  pm  thick)  on  a  c-plane  sapphire  substrate. 
Solid  line  represents  the  isotropic  sapphire  or  GaN  layer;  dash  line  represents  the 
anisotropic  layers. 

Fig.  2.  A  simulation  sample  structure  with  the  surface  over-layer  for  VASE  data 
generation  the  surface  roughness  is  described  by  a  EMA  layer  with  50%  voids. 

Fig.  3.  Surface  roughness  effect  observation  through  VASE  data  simulations  (angle 
dependence).  Generated  VASE  data  of  Fig.  2  structure  at  2  eV  and  5  eV  by  isotropic 
model  with  different  over-layer  thickness,  (a)  'f  as  a  function  of  angle  of  incidence  at  5 

eV;  (b)  A  as  a  function  of  angle  of  incidence  at  5  eV;  (c)  VF  as  a  function  of  angle  of 
incidence  at  2  eV;  (d)  A  as  a  function  of  angle  of  incidence  at  2  eV. 

Fig  4.  Surface  roughness  effect  observation  through  VASE  data  simulations  (energy 
dependence).  Generated  VASE  data  as  a  function  of  photon  energy  at  60°  and  80°  angles 
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of  incidence  by  isotropic  model  with  different  over-layer  thickness,  (a)  'F  spectra  with  0 


and  4  nm  thick  over-layer;  (b)  A  spectra  with  0  and  4  nm  thick  over-layer. 

Fig.  5.  Surface  roughness  effect  observation  through  reflectance  data  simulations  (both 
angle  and  energy  dependence),  (a)  Generated  reflectance  data  as  a  function  of  energy  at 
60°  and  80°  angles  of  incidence;  (b)  Generated  reflectance  data  as  a  function  of  angles  of 
incidence  at  2  eV  photon  energy  (below  band  gap);  (c)  Generated  reflectance  data  as  a 
function  of  angles  of  incidence  at  5  eV  photon  energy  (above  band  gap). 

Fig.6.  The  three-phase  model  of  the  MBE  grown  GaN  on  c-plane  sapphire  sample. 

Fig.  7.  Experimental  VASE  data  and  the  best  model  fit  of  a  GaN/c-sapphire  sample  at 
three  angles  of  incidence.  The  circles  and  solid  line  represent  the  experimental  data  and 
the  best  fit,  respectively. 

Fig.  8  Ordinary  and  extraordinary  dielectric  functions  determined  from  the  VASE  data 
analysis.  The  solid  line  represents  the  ordinary  dielectric  functions  (ex).  The  dash  line 

represents  the  extraordinary  dielectric  functions  (ey). 
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Fig.  9.  Comparison  of  the  experiment  results  and  the  calculated  £1  value  for  both  the 

ordinary  and  extraordinary  dielectric  functions  of  GaN.  The  solid  dots  and  circles 
represent  experiment  results,  while  the  lines  represent  calculated  values  by  K-K  fit. 


17 


Table  1 .  The  best  fitting  parameters  for  GaN  extraordinary  dielectric  functions. 


CP  contribution 

Aj  (Ajx) 

Ej  (EjX) 

Ej  (rjx) 

Pjx 

e<°> 

38.5  (eV)15 

3.44  (eV) 

0.27  (eV) 

E(°x) 

1.32  (eV) 

3.42  (eV) 

0.18  (eV) 

-1.55 

e(l) 

2.67 

6.96  (eV) 

0.003  (eV) 

£(lx) 

4.77  (eV) 

6.94  (eV) 

1.38  (eV) 

-0.01 

0.24 
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Table  2.  Kramers-Kronig  fitting  parameters  for  GaN  ordinary  dielectric  functions. 


Fitting  Parameters 

Ei 

Aj 

0  offset 
£i 

Pole  #1 

7.491 

2.001 

1.04 

Pole  #2 

14.092 

24.6 
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Table  3.  Kramers-Kronig  fitting  parameters  for  GaN  extraordinary  dielectric  functions. 


Fitting  Parameters 

Ei 

Aj 

0  offset 
£i 

Pole  #1 

7.4883 

3.805 

1.05 

Pole  #2 

13.534 

21.123 
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